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We consider instability of the Friedmann world model to the second-order in perturbations. We 
present the perturbed set of equations up to the second-order in the Friedmann background world 
model with general spatial curvature and the cosmological constant. We consider systems with the 
completely general imperfect fluids, the minimally coupled scalar fields, the electro-magnetic field, 
and the generalized gravity theories. We also present the case of null geodesic equations, and the one 
based on the relativistic Boltzmann equation. In due stage a decomposition is made for the scalar-, 
vector- and tensor-type perturbations which couple each other to the second-order. Gauge issue is 
resolved to each order. The basic equations are presented without imposing any gauge condition, 
thus in a gauge-ready form so that we can use the full advantage of having the gauge freedom 
in analysing the problems. As an application we show that to the second-order in perturbation 
the relativistic pressureless ideal fluid of the scalar-type reproduces exactly the known Newtonian 
result. As another application we rederive the large-scale conserved quantities (of the pure scalar- 
and tensor-perturbations) to the second order, first shown by Salopek and Bond, now from the exact 
equations. Several other applications are made as well. 

PACS numbers: 04.50. +h, 04.62. +v, 98.80.-k, 98.80.Hw 



11 the linear perturbations are so surprisingly simple that a perturbation analysis accurate to second order 
may be feasible ..." 

Sachs and Wolfe (1967) 
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I. INTRODUCTION 

We consider instabilities of the spatially homogeneous and isotropic comological spacetime to the second-order in 
perturbations. The relativistic cosmological perturbation plays fundamental roles in the modern theory of large- 
scale cosmic structure formation. The original analysis of linear perturbations based on Einstein gravity with a 
hydrodynamic fluid was made by Lifshitz in 1946 [1] in an almost complete form. Due to the extremely low level 
anisotropics of the cosmic microwave background (CMB) radiation, the cosmological dynamics of the structures in 
the large-scale and in the early universe are generally believed to be small deviations from the homogeneous and 
isotropic background Fricdmann world model [2]. The conventional relativisitic cosmological perturbation analysis 
considers such deviations small enough so that one could treat them as linear. The linear perturbation theory works 
as the basic framework in handling the cosmological structure formation processes. Recent observations of the CMB 
anisotropies in all sky by the WMAP satellite and others [3], for example, assure the validity of the basic assumptions 
used in the cosmological perturbation theory, i.e., the linearity of the relevant cosmic structures [4]. 

Still, as the observed relatively small-scale structures are apparently nonlinear, the gravitational instability based 
on the pure linear theory is not enough for the complete picture. It is agreed that such small-scale nonlinear structures 
could be handled by Newtonian gravity often based on numerical simulations. The current paradigm of large-scale 
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structure generation and evolution processes is based on an underlying assumption that linear processes dominated 
until nonlinear processes take over on subhorizon scales in the Newtonian regime. Thus, it seems this paradigm of our 
understanding of the origin and evolution of the large-scale structures is rather satisfactory within the linear theory 
concerning the regimes where the relativistic gravity theory is needed. It is well known that in the linear theory 
there can be no structure formation. In fact, this "no structure formation" in the scenario is precisely why we were 
successful in describing the structure generation and evolution processes in a simple manner as we describe below. In 
the standard scenario, the initial conditions (seed fluctuations) generated from quantum fluctuations are imprinted 
into ripples in the spacetime, and its spatial structures are preserved as the raw large-scale structures (i.e. are as 
yet unaffected by nonlinear processes). This is a trait which can be traced to the linearity assumption we adopt. It 
gives a simple, but fictitious system. We could, perhaps, fairly describe the current situation as: the linear paradigm 
is not inconsistent with observations, especially with the low level of observed anisotropy in the CMB [3]. However, 
we should remember that the actual equations we are dealing with, both in the gravity and the quantum, are highly 
nonlinear. It forms an intrinsically complex system. 

We can decompose the perturbations into three different types: the scalar-type (associated with the density con- 
densation), the vector- type (rotation) and the tensor- type (gravitational wave) perturbations. To the linear-order 
perturbation, due to the high symmetry of the background space, these three types of perturbations decouple from 
each other and evolve independently. Both (the linearity and the homogeneous-isotropic background) conditions 
are necessary to have natural descriptions of the three types of perturbations independently. We will see how cou- 
plings occur to the second-order perturbations; for the couplings in a simplest spatially homogeneous but anisotropic 
spacetime, see [5]. 

Another aspect of the simple nature of the linear processes in relativistic gravity theories is characterized by 
conservations (in expanding phase) of certain amplitudes in the super-horizon scale where we naively anticipate the 
independent evolution of causally disconnected regions. On superhorizon scales this conserved character is presented 
by an equation of the form 

*(x,t)=C(x), (1) 

which applies for both the scalar- and tensor-type perturbations; for the vector-type perturbation the angular momen- 
tum is conserved in all scales. C(x) is an integration constant of the integral form general solutions available in the 
large-scale limit, see eqs. (311,332,333). The coefficient C(x) contains information about spatial structure which will 
eventually grow to the large-scale structure and the gravitational wave background. It can be considered as an initial 
condition for each perturbation variable which is preserved during the linear evolution. Whether a similar conserved 
variable can exist even in a nonlinear analysis is apparently an interesting question: for the presence of such variables 
to the second-order in perturbation see §VIID-VIIF. In the analyses of the large-scale structures in the linear stage, 
the simple behavior of the conserved variables is practically important. In fact, if we know C the behavior of all the 
other variables can be determined through linear algebra. Using the conserved quantity one can trivially relate the 
currently observable (or deducible) linear structure directly to the initial state of the structure at the early universe; 
probably, just after the scale effectively becomes the large-scale during the hypothetical early acceleration (inflation) 
stage. Of course, the underlying assumption for all of these results is the applicability of the linear analysis. As long 
as this assumption is valid, the initial condition is imprinted onto the large-scale structure and is preserved until the 
nonlinear effects become important. Although this is a big advantage, in a sense this is very consistent with the fact 
that no structure formation occurs in the linearized system. 

The linear perturbation theory is currently well developed; see [6-11]. Although the observations do not desparately 
demand to go beyond the linear theory, the second-order perturbation theory is a natural next step in the theoretical 
investigations. The second-order perturbation theory, if well developed, will have important implications on our 
understanding of the large-scale structure formation processes. Not only the structures we discover are nonlinear, 
according to the gravitational instability there should occur (perhaps smooth) transitions from the linear to the 
nonlinear ones. Even in the theoretical point of view, in order to know the limit of the linear perturbation theory 
we need the behavior of perturbations beyond the linear theory. It is not possible to know the limit of linear theory 
within the context of the linear theory. It is yet unclear whether the second-order 'perturbation theory' will provide 
an answer to such a question, but we expect it could provide a better perspective on the problem than the simple 
linear theory. There will be more practical applications as well, like investigating the non-Gaussian signature in the 
inflation generated seed fluctuations which could have left a detectable signature in the CMB anisotropies and the 
large-scale structures. Other possible situations where the trans-linear analyses might be useful are summarized in 

§vm. 

Now, we discuss the gauge issue present in the relativistic perturbation theory briefly. Since the unperturbed 
background spacetime is spatially homogeneous and isotropic, to the linear-order the ambiguity caused by spatial 
gauge (coordinate transformation) freedom does not play a role [7]. Thus, to the linear-order it is appropriate to 
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write the perturbed set of equations in terms of natural combinations of variables which are invariant under spatial 
gauge transformation. However, the temporal gauge freedom can be conveniently used in analyzing various aspects of 
perturbation problem in the Friedmann background. There are infinitely many different ways of taking the temporal 
gauge (choosing the spatial hypersurface) conditions, and we can identify several fundamental gauge choices [7,12]. 
Except for a widely used temporal synchronous gauge fixing condition (Sgoo = 0), each of the other gauge conditions 
fixes completely the temporal gauge freedom. Thus each having its own corresponding gauge invariant formulation. 
In our study of the linear theory we found that some particular gauge invariant variables show the correct Newtonian 
behaviors. A perturbed density variable in the comoving gauge and a perturbed potential variable and a perturbed 
velocity variable in the zero-shear gauge most closely resemble the behaviors of the corresponding Newtonian variables 
[13]. Also, the scalar field perturbation in the uniform-curvature gauge most closely resembles the scalar field equation 
in the quantum field in curved space [14,15]. Since the gauge conditions mentioned allow no room for remaining gauge 
mode the variable in a given such gauges has the unique corresponding gauge invariant combination of variables. Thus, 
the variable in such a gauge is equivalent to the corresponding gauge invariant combination. 

In the gauge theory it is well known that proper choice of the gauge condition is often necessary for proper handling 
of the problem. Either by fixing certain gauge conditions or by choosing certain gauge-invariant combinations in the 
early calculation stage we are likely to lose possible advantages available in the other gauge conditions. In order 
to use the various gauge conditions as advantage in handling cosmological perturbations we have proposed a gauge- 
ready method which allows the flexible use of the various fundamental gauge conditions. The strategy is that, in 
order to use the various available temporal gauge conditions as advantage, we had better present the basic equations 
without choosing any temporal gauge condition, and arrange the equations so that we could choose the fundamental 
gauge conditions conveniently. In this work we will elaborate further the gauge-ready approach to the second-order 
perturbations. Our gauge-ready strategy, together with our notation for indicating the gauge-invariant combinations, 
allows us to use the gauge freedom as the advantage in analysing various given problems. We follow the wisdom 
suggested by Bardeen in 1988 "the moral is that one should work in the gauge that is mathematically most convenient 
for the problem at hand" [7] . 

As long as we are taking perturbative approach the gauge issue in the higher-orders can be resolved similarly as 
in the linear theory. To the second-order we will identify two variables which can be used to fix the spatial gauge 
freedom. One gauge condition completely removes the spatial gauge mode, whereas the other condition does not; i.e., 
in the latter case even after imposing the gauge condition there still remains a degree of freedom which is a gauge 
mode (a coordinate artifact). We call this incomplete gauge condition the -B-gauge (goa = where a is a spatial 
index) whereas the other complete condition is called the C-gauge. To the second-order we can identify the same 
several temporal gauge fixing conditions. Again, except for the synchronous gauge each of the other gauge conditions 
removes completely the temporal gauge modes. 

It is amusing to note that the classic study by Lifshitz [1] adopted the synchronous gauge condition which is a 
combination of the temporal synchronous gauge condition and the spatial .B-gauge condition, thus failing to fix both 
the temporal and the spatial gauge modes completely. This has caused some prevalent errors in the literature based 
on the synchronous gauge: see the Appendix of [13]. However, we note that these errors are simple algebraic ones 
probably caused by slightly more complicated algebra due to the presence of the gauge mode after the synchronous 
gauge fixing. We would like to emphasize that the gauge condition should be appropriately used according to the 
character of each problem at hand. We have such a freedom because Einstein's gravity theory might be regarded 
as a gauge theory [16]. In this sense, although the temporal synchronous and the spatial i3-gauge conditions do not 
remove completely the gauge modes, often even these conditions could possibly turn out to be convenient in certain 
problems. Since physically measurable quantities should be gauge invariant we propose to use the gauge conditions 
in this pragmatic sense. 

In a classic study of CMB anisotropy in 1967 Sachs and Wolfe have mentioned that "the linear perturbations are 
so surprisingly simple that a perturbation analysis accurate to second order may be feasible ..." [17]. In this work, 
we will present the basic formulation of the second-order perturbation of the Friedmann world model in details. We 
will present the basic equations needed to investigate the second-order perturbation in a rather general context. We 
will consider the most general Friedmann background with K and A. We will consider the most general imperfect 
fluid situation. This includes multiple imperfect fluids with general interactions among them. We will also include 
minimally coupled scalar fields, a class of generalized gravity theories, the electromagnetic fields, the null geodesic, and 
the relativistic Boltzmann equation. In order to use the gauge fixing conditions optimally we will present the complete 
sets of perturbed equations in the gauge-ready form. In this manner, as in the linear theory, we could easily apply the 
equations to any gauge conditions which make the mathematical analyses of given problems simplest. Our formulation 
will be suitable to handle the nonlinear evolutions in the perturbative manner. Ours will be a useful complement to 
the other methods suggested in the literature to investigate the trans-linear regimes. In the Discussions we summarize 
the related studies, the different methods we could employ for further applications, and the cosmological situations 
where our formultion could be applied fruitfully. Although we will present some trivial applications in the later part, 
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the main applications are left for future studies. 

In §11 we summarize the basic equations of the Einstein's gravity expressed in the ADM (3 + 1) formulation and in 
the covariant (1 + 3) formulation. In §111 we introduce our definition of the metric and the energy-momentum tensor 
to the second-order in perturbation and present some useful quantities appearing in the ADM and the covariant 
formulations. In §IV we present the complete sets of perturbed equations up to the second-order in the Friedmann 
world model. We consider the general spatial curvature and the cosmological constant in the background. We consider 
systems with the completely general imperfect fluids. Such a general formulation can be reinterpreted to include the 
cases of the minimally coupled scalar fields, the electro-magnetic field, and even the generalized gravity theories. We 
present the complete sets of equations for these systems as well. We also present the case of null geodesic equations, 
and the one based on the relativistic Boltzmann equation. In §V we introduce decomposition of perturbations to 
three different types and show how these couple to each other to the second-order. All equations up to this point 
are presented without introducing any gauge condition. Thus, the equations are presented in the most general forms, 
and any suitable gauge conditions can be easily deployed to these equations. In this sense, our set of equations is 
in a gauge-ready form. In §VI we address the gauge issue, and show that the gauge issue can be resolved to each 
perturbation order, just like the case in the linear perturbation. We implement our gauge-ready strategy to the 
second-order in perturbations. In §VII we make several applications. In §VIII we summarize the main results and 
outline future applications of our work. 

As a unit we set c = 1. 



II. BASIC EQUATIONS 
A. ADM (3 + 1) equations 

The ADM (Arnowitt-Deser-Misner) equations [18] are based on splitting the spacetime into the spatial and the 
temporal parts using a normal vector field n a . The metric is written as (we put tilde on the covariant variables): 

goo = -A 2 + N a N a , g 0a = N a , g a p = h a p, 

g oa = -N- 2 , g aa =N- 2 N a , g af3 = h af3 - N~ 2 N a N 13 ' , (2) 

where N a is based on h a p as the metric and h af3 is an inverse metric of h a p- Indices a,b, . . . indicate the spacetime 
indices, and a, /3, . . . indicate the spatial ones. The normal vector h a is introduced as: 

n ee —N, h a ee 0, n° = N' 1 , h a = -N^N". (3) 

The fluid quantities are defined as: 

E ee h a hi,T ab , J a = — n&T^, S a p ee T a p, S ee h af3 S a p, S a p ee S a p — -h a pS, (4) 

where J a and S a p are based on h a p. The extrinsic curvature is introduced as 

K a p ee (JV q:j8 + Np :a - h a p fi ) , K=h af3 K a p, K a[j = K al3 -^h al3 K, (5) 

where K a p is based on h a p. A colon ':' denotes a covariant derivative based on h a p. The connections become: 

f So = ^ {No + A Q A Q - K af3 N a NP) , f ° Qa = 1 (N a - K af3 N?) , T% = ~K aP , 

fg = ^ a (-N - NjiN t3 + KfjjN^N 7 ) + NN' a + N a - 2NK af3 N p + N a:(3 N fil 

% = -^N,pN«-NK% + N a ^ + ^N<*WK Pl , f% = T {h) ^ + ^N a K^ (6) 

where T^jf is the connection based on h a p as the metric, r^^j ee \h aS (hps,-y + hs^^p — h^j). The intrinsic 
curvatures are based on h a p as the metric: 

T>(h)a _ -p(h)a _ p(h)a . -p(h)e -p(h)a _ -p(h)e -p(h)a 
11 PjS — 1 06,-y 1 /37,5 + i pS L 7« 07 5e> 

Ri h h R{h) \^ R^^h^R^, R%^R { S~\h a pR^. (7) 
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A complete set of the ADM equations is the following [6] . 
Energy constraint equation: 



R (h) = K a(3 K aP - \k 2 + WttGE + 2A, (8) 



where A is the cosmological constant. 
Momentum constraint equation: 



Trace of ADM propagation equation: 



K* :f) - \K. a = 8nGJ a . (9) 



A qN- 1 - K^ a N a N~ 1 + N :a a N^ 1 - K afi K a[j - \k 2 - 4nG (E + S) + A = 0. (10) 
Tracefree ADM propagation equation: 

= KKp - (N :a - ^N : \j N- 1 + R [h)a fj - 8ttGS%. (11) 



Energy conservation equation: 



oA" 1 - E^ a N a N~ 1 -k(e+ ^S^j - S afi K afi + N~ 2 (N 2 J a ) -q = 0. 



(12) 



Momentum conservation equation: 

J^oN- 1 - J^N^N- 1 - JpNP, a N- x - KJ a + EN^N- 1 + Sf :/3 + S^N^N -1 = 0. (13) 

B. Covariant (1 + 3) equations 

The covariant formulation of Einstein gravity was investigated in [19,20]. The 1 + 3 covariant decomposition is 
based on the time-like normalized {u a u a = —1) four-vector field u a introduced in all spacetime points. The expansion 
(9), the acceleration (a a ), the rotation (ui a b), an d the shear (<r a b) ar e kinematic quantities of the projected covariant 
derivative of flow vector u a introduced as 

h c a hfu c . d = h c [a h b] ii c . d + h c (a h^u cid = Lj ab + 6 ab = u a , b + a a u h , 

Oab = d ab - ^Ohab, = u a . a , a a = u a = u a . b u b ', (14) 

where h a b = g a b + u a u is the projection tensor with h a bU = and h a a = 3. An overdot with tilde ' indicates a 
covariant derivative along u a . We have 

U a ;b = Uab + &ab + ^6h ab - CL a U b . (15) 

We introduce 

,~, a — ^ a bcd~ ~ ~ ~ r, c ?.d ~2 — ^ ~,ab ~ r, a ,~. ~2 — ^ ~ab ~ ! ~\ a\ 

U — u b^cd, Uab = VabcdU u , UJ = -UJ UJ ab — UJ Ul a , a = -a a ab , (lb) 

where uj a is a vorticity vector which has the same information as the vorticity tensor uj ab . We have fj abcd = ff abcd ^ with 
^1234 _ i/yf^g- indices surrounded by () and [] are the symmetrization and anti-symmetrization symbols, respectively. 
Our convention of the Riemann curvature and Einstein's equation are: 

Ua-bc - Ua;cb = U d R d abc , (17) 

Rab - \R~gab = SnGfab - Ag ab . (18) 
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The Weyl (conformal) curvature is introduced as 

1 / \ R 

Cabcd = Rabcd — \9acRb d + 9bdRac — 9bcRad - 9adRbc) + (9ac9bd ~ 9ad9bc) ■ (19) 

The electric and magnetic parts of the Weyl curvature are introduced as: 

E ab = C acbd u c u d , H ab = ^f) ac ef C ef b d u c u d . (20) 

The energy-momentum tensor is decomposed into fluid quantities based on the four-vector field u a as 

f ab = jlu a Ub + p {g a b + u a u b ) + q a u b + q b u a + fr ab , (21) 

where 

u a q a = = u a n ab , iab = ^ba, K = °- ( 22 ) 

The variables jl, p, q a and ir ab are the energy density, the isotropic pressure (including the entropic one), the energy 
flux and the anisotropic pressure based on u a -frame, respectively. We have 

fi = f ab u a u\ p=lf ab h ab , q a ee -f cd u c h d a , n ab = f cd hfh d -ph ab . (23) 

The specific entropy S can be introduced by TdS = de + prdv where e is specific internal energy density with 
/i = g(l + e), Pt the thermodynamic pressure, v = 1/g the specific volume, and T the temperature. We have the 
isotropic pressure p = pr + tt where n is the entropic pressure. Using eq. (26) below we can show 

gTS = - (#0 + n ab a ab + q a -, a + Q a a a ) • (24) 



entropy flow density [19] we can derive 



Thus, we notice that tt, w and q a generate the entropy. Using a four-vector S a = gu a S + ~q a which is termed the 



S a . a = ~{f >a + fa a ) q a - I (tt6 + 7r ab a ab ). (25) 

The covariant formulation provides a useful complement to the ADM formulation. We summarize the covariant 
(1 + 3) set of equations in the following. For details, see [19,20] and the Appendix in [21]. 
The energy and the momentum conservation equations: 

/*+(£+ P) & + ^ h Oab + q a :a + Q^a = 0, (26) 

(/i + p) ~a a + h b a (p tb + 7Tg. c + + (ti ab + o ab + Uh ab ^j q b = 0. (27) 



Raychaudhuri equation: 



Vorticity propagation: 



9 + \o 2 - d a . a + 2 (a 2 - u 2 ) + 4ttG (p + 3p)-A = 0. (28) 



h a b tb b + ^9Cj a = d a b u b + X -ff bcd u b d c , d . (29) 



Shear propagation: 

- ~ /~ \ 2 ~ 1 ~ 

K.K [Zed - a( c; d)J - d a d b + Lu a Lu b + a ac a b + -8a ab - -h ab (a) 2 + 2a 2 - d c . c ) + E ab - AitGTi ab = 0. (30) 

Three constraint equations: 
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Kb (V c . c - a bc . c + p b ^j + ((b ab + a ab ) ~a b = 8nGq a , (31) 

u a , a = 2£j b ~a b , (32) 

H ab = 2~a (a u b) - h c J4 (cD (c e;/ + a (c e;/ ) r,d) gef u 9 ■ (33) 

Four quasi-Maxwellian equations: 

hfh c d E bd . c - f, abcd u b alH de + SH^o = AnG (f^fti - ;c - 3a>°i/ + a a b q b + ^a b - \k a ^j , (34) 

h a b h c d H bd , c + fj abcd u b a e c E de - 3E^ b = 4ttG{2(p + p) Cu a + fj abcd u b [q c . d + n ce (Cj e d + a e d )]} , (35) 
h a c h b d E cd + (Hi e h { f - Ta d Hi a ^ f] b ^ de u c + h ab b cd E cd + 6E ab - M c a (3^ c + £ b ' c ) = 4ttG[ - (fi + p) o ab 

-2a^ 6) - hi a h b J (q c - d + ^) - (w c (a + ^) ^ ~ \^ ab + \ {f , c + a c q c + n cd a cd ) h ab \ , (36) 
h a c h b d H cd - (ElJi^ - 2a d E^ f] b)cde u c + h ab d cd H cd + 6H ab - H^ a (zd b)c + £ fc)c ) 

(q e d { d a - ^ f d . e h [ f) V b)cde u c + h ab Cj c q c - 3a}(°g 6 )] . (37) 
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Evaluated in the normal-frame eqs. (26,27,28,30,31) reproduce eqs. (12,13,10,11,9) in the ADM formulation. 

Now, we take the normal-frame vector, thus u a = n a with n a = 0, thus £j ab = 0. The trace and tracefree parts of 
the Gauss equation give [21]: 

R (3) = 2 (-^d 2 +° 2 +8nGfi + A^j , (38) 

R-ab ~ ^^Kb = h c J4 (-Zed - 0v cd + a (c;d) ) + a a a b - ^h ab a c . c + 8irGn ab , (39) 

where R a 3 b ] and R^l are the Ricci and scalar curvatures of the hypersurface normal to h a ; for an arbitrary vector V a 
we have 



5(3) frb _ 
H abcd V 



2v[fvyy a E2^^»v [e (^ ] / i ;v A y 1 ), r%> = h cd R% db , RW = h ab R%. (40) 

From this we have 

Raid = KHkh h d Refgh ~ B ca 6 db + hAd, (41) 

which is the Gauss equation. We can show that R^p^s = ^afiyS- Equation (39) follows from cq. (30) evaluated in the 
normal-frame. Using eqs. (19,20) we can show that eq. (39) reproduces eq. (11) in the ADM formulation. Equation 
(38) gives cq. (8) in the ADM formulation. 

Compared with the ADM equations in (8-13) part of the covariant equations in (29,32,33,34-37) look new. In the 
normal-frame eqs. (29,32) are identically satisfied; using eqs. (14,6,3) we can show 

5 a = (InJV) ia) (42) 

thus «[/3 ;7 ] = 0. Still, eqs. (8-13) provide a complete set. These additional equations in the covariant forms should 
be regarded as complementary equations which could possibly show certain aspects of the system better. In our 
perturbation analyses we will use parts of these equations as complement ones. Although the covariant set of equations 
is based on the general frame vector, this does not add any new physics which is not covered by the normal-frame 
taken in the ADM formulation, sec §111 E. 

The covariant equations for the scalar fields, generalized gravity, electromagnetic field, null geodesic, and Boltzmann 
equation will be introduced individually in the corresponding sections later. 

C. Multi-component situation 

In the multi-component situation we have 



8 



fab = T(l)*b; f (l) l b = J (i)o , Y h)° = °' ( 43 ) 
I I 

Based on the normal-frame vector, we have 

A = X^(')> P = X^(0' Qa = ^2q(i) a , ir ab = ^2n (l)ab . (44) 
ill i 

The ADM formulation is based on the normal- frame vector u a = h a . The ADM fluid quantities in eq. (4) correspond 
to the fluid quantities based on the normal-frame vector as 

E = ji, S = 3p, J a =q a , S a = 7t a0 . (45) 

From eq. (4) or eq. (45) we have 

i « ; j 

Equation (43) gives 

i^oiV- 1 - E^N'N- 1 - K (e^ + - Sg>tf a/J + AT 2 (a^) ^ - (/(i)o - ^)a^ a ) , (47) 

JWcoN- 1 - Ja^pN^N- 1 - J^NP.^N- 1 - KJ (i)a + E^NaN- 1 + 5 (i) f :/3 + JV^JV" 1 = J (i)a . (48) 

The ADM equations in eqs. (8-13) remain valid, with the above additional equations of motion for the individual 
component. Thus, in the multi-component situation eqs. (8-13,46-48) provide a complete set. 

III. PERTURBED QUANTITIES 
A. Metric and connections 

We use the following convention for the metric variables: 

goo = -a 2 (1 + 2A) , g 0a = -a 2 B a , ~g a p = a 2 {g^j + 2C Q/3 ) , (49) 

(31 

where A, B a and C a p are perturbed order variables and are assumed to be based on g y a ^ as the metric. To the 
second-order, we can write the perturbation variables explicitly as: 

A = AM+AW, Caf, = C<$+C$. (50) 

As we are interested in the perturbation to the second-order, as our ansatz, we include up to second-order (quadratic) 
terms in the deviation from the Friedmann background. This can be extended to any higher-order perturbation 
as long as we take the perturbative approach where the lower-order solutions drive (work as sources for) the next 
higher-order variables. Thus, in this work we ignore the terms higher than quadratic (second-order) combination of 
the perturbed metric (A, B a , C a p), the perturbed fluid quantities (Sfj,, Sp, Q a , H a p) to be introduced in eq. (72), 
the perturbed field (5(j)) to be introduced in eq. (Ill), etc. 

The inverse metric expanded to the second-order in perturbation variables is: 

9°° = \(-l + 2A- 4A 2 + B a B a ) , 
g 0a = 4 (-B a + 2AB a + 2B C al3 ) , 

~ g a P = -L ( 5 (3)«/3 _ 2C af3 - B a B f3 + AC*C ri ~<^j . (51) 

The connections are: 
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fgo = - + A' - 2AA' - A^B a + B a (B 



a, , ^ B c 

a 



T a = A 1 " - B Q ' - — B a + A'B a - 2A^C afi + 2C% { B f3/ + —B 



1 a/3 



1 0/3 



ffi? - + B (aW + C' af} + 2^-C af3 



.(3) 



+4S ( 4 ^ 2 - B , B1 ) - 2A ( B WV + C" a0 + B, (2C^ - C a f) , 



fia _ p(3)a I a „( 3 ) I oAia I" n^a ( cy^S f~i \S\ 

l (3~<- 1 fi~, + —9p~t B + 2C (/3|7) " C /37 ~ 2C <5 V 2C (/3|7) ~ C /3 7 J 



(52) 



where a vertical bar indicates a covariant derivative based on g y a ^ . An index indicates the conformal time 77, and a 
prime indicates a time derivative with respect to rj. The components of the frame four-vector u a are introduced as: 



u=-ll — A + - A + -V a V a - V a B a , u a = -U Q , 
a \ 2 2 / a 

wo = -a + A - X -A 2 + \v a v}j , w a = a (U Q - B a + AB a + 2V (i C a(i ) , 



(53) 



where V a is based on . 

B. Normal-frame quantities 

The normal- frame vector h a has a property h a = 0. Thus we have 



i° = - (l- A+-A 2 - -B a B a ) , n Q ee - (B a - AB a - 2B fi CS) 



ih, = '< I 1 + A - ^A 2 + \B a B a ) »„ = 0. 



(54) 



Using eqs. (2,3) the ADM metric variables become: 

a = -a 2 B a , N a = — B a + 2B fj C%, 
K fi = a 2 (g% + 2C a ,) , = 1 (gW"P 2C^ + 1C«C^) . 



V = „ I • .1 \a 2 • -/> "■■) , 



The connection becomes 



<hh _ F (3)7 



a/3 



a/3 



The extrinsic curvature in eq. (5) gives: 



K, 



a/3 



K = — 
a 



+ *<«|0 + C' ap + 2 a -C a ^ (1 - A) + \ a -g§ (3A 2 - B,Bi) B, ( 2 ^ |/3) - cj") 



(55) 



(56) 
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K afj = -a 
1 



| (B {am + C' ap ) {I -A)- B, (2Cj am C'Jt) \c afj (ST |7 + C]' 



(3) 



The intrinsic curvature in eq. (7) becomes: 



(57) 



R (h ) = 



a/3 + (^g i ' 3 ^ S — 2C 7(5 ^ (Cs a \i3 7 + C$p\ al — C Q/3 |,5 7 — Cs^afj) + 2C' jS ^C 7 s\a 
~ ^ °5|7 ~~ U 1\SJ [y<x\l3 + C P\a - ^afj J ~ {^ a \ 7 + Sl« ~ °"7 J \ C /3\S + U 6\/3 ~ °/3<5 J' 



where 



/3 7 5 g 

with a normalized K(= 0, ±1), the sign of the background three-space curvature. Thus, 

f>{h)a _ 1 J n a~t 1 /Yll" f~ia\l /~a\a ^ (3) 

^~^21 107 ^ 7 ~ " 7 7 ^"3^ C /3 



(58) 



(59) 



^°5|/3 7 + ^5/3 7 - C /3|5 7 -°5 7 ^ 7 |/3« + C /3| 7 5 ~ °/37 <5 ~ °« 17/9 J + 3 K °7 °, 

- ( 2 ^l7 - ^1 *) + C T ~ ^) + ^l/^ 1 * + 2C"^ (C^l, - C^| 7 ) 

\r&g< + 2cf s 2cf s + ac^ s (-c; ]Se c 7|e4 + c l5 |e e + c^ s ) 

+±RV>C*(P S - (2Cl le ~ Ct\ S ) (2C 5 h - C^ s ) + (3C 75 , e - 2C lA& 
It is convenient to have 



3^ 



B a ^=B 



1/37 ~ " 1 7/3 



R {3)a xa B 5 , B, 



S0j° > ^a|/3 7 - ^a| 7 /3 



(3)5 



a/37 



(60) 



(61) 



C. General ({t a )-frame quantities 

To the second-order perturbation, using cqs. (49,51,53), the kinematic quantities in eq. (14) become: 
Kp = a 2 \gW + 2C aP + (V a - B a )(V p - Bp)] , h 0a - -a 2 (V a + AV a + 2C al3 V f} ) , V> - a 2 V a V ai 



~e = a - 1 



3— (1 - A) + V a \ a + C% + — A 2 + B a B' a - 2C afi C' afj - AC% 



2 a 



~V a 'B a + V- (v> + 3 -^- Va -B' a - Z a -B a + A, a + Cf |a ) 



a a = A, a + V^-B' a + - (V a - B a ) + A'B a + A -2A, a + 2B' a + 2-B a -V' a - -V a 
a V a a 



(62) 



(63) 
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+V 3 (V a \ fi + B fM - B ay} ) + 2C a(i (vf» + ^V) + 2C' a pVP, ~a a = -V a a a , 
-V [a V/ 3] + B [a Vh + V [a B' 0] + 2V 1 C^ am + 2C l[a V m ) , Cu 0a = V^u, a p, Co oa = 0, 



a a p = a 



/3) 



+V (a Aj3) + V 7 C a ^ - AC' a!i + 2C 5{a V 5 m - -C a p (V~< h + C*) 

(vW] - VB'^ - V^'B 1 + B~>B' y + VM, 7 + V S C^ S - AC]' - 2C^C'^ 

&0a = -V^O-ap, (Too = 0. 

In the normal-frame we have u a = 0, thus V a = B a — AB a — 2B^C a p. In this frame we have 

h a p = a 2 (g^j + 2C af ^j , h 0a = -a 2 B a , h oa = a 2 B a B a , 
6 = -K, 

~a a = (\nN) a = ^A-A 2 + ^Bf } Bp^j , a = -B a A, a , 
& a f3 = —K a p, a 0a = B^K a p, CT o = 0, 

Uab = 0. 



(64) 



(65) 



(66) 

(67) 
(68) 
(69) 

(70) 
(71) 



In this frame we have 9 = —K and a a p — —K a p. These are natural because K and K a p are the same as negatives 
of the expansion scalar and the shear, respectively, of the normal-frame vector field. 



D. Fluid quantities 

To the perturbed order we decompose the fluid quantities as: 

p, = fj, + Sfi, p = p + Sp, q a =aQ a , n a p = a 2 H a /3, 
,(3) 



(72) 



where Q a and H a p are based on g a p . In the Friedmann world model we have ft = n and p = p and zeros for the other 
fluid quantities. We have 



n a _ 2C af3 n a p = 0, 



(73) 



which follows from 7r^ = or S™ — 0. The perturbed order fluid quantities 5[i, Sp, Q a and H a p in eq. (72) can be 
expanded similarly as in equation (50) 

6» = 6^+6^ 2 \ 8p = 5pW+5 P <», Q a = QP + Q™, n Q/3 = n$ + ng> . (74) 

In the multi-component situation, from eqs. (44,72) we set 

i i 

^ = 5Z (5/i (0' s P = ^2 s P(i), Q a = ^2Q(i) a , n Q/3 = ^2u (l)a p. (75) 
iii i 

Thus, from eq. (45) the ADM fluid variables become: 
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E = fj, + (5/i, 
S = 3(p + 6p), 

J a = aQ a , J a = 1 (Q a - 2C af3 Qp) , 

s aP = a 2 n a p, s% = n<$-2C a m Pj , s af3 = ± (n a(} -4C^ Q n«) . (76) 

From eqs. (22,72) we have 

q a = aQ a , <jo = —aQ a B a , 

TT a p = a 2 U af3 , TT 0a = —a^HapB 13 , 7foo = 0. (77) 

From cqs. (21,72,54) we have 

Too = a 2 [fi + Sfi + 2(j,A + 2SfiA + (/i + p) B a B a + 2Q a B a ] , 
{Qa + pB a + 5 P B a + AQ a + n Q/3 B /3 ) , 

f a p = a 2 (pg™ + dpgfl + U a p + 2pC a p + 25pC a p) . (78) 

From eqs. (4,76) we have 

f i + 5n = f ab n a n b , p + 5p=\h al3 f a p, Q a = --f ab n\ U afj = \ (f a p - h a[i p) . (79) 

O CI CL \ / 

Equation (24) gives 

gf§ = ttK+ \n a PR a p - 2-Q a A. a - -Q a , + 2- (C^Qp), - -C^Q fj . (80) 
a 2 a a 1 a '\ a a lP 

For the interaction terms in eq. (43) we set 

l(i)0 = l(i)0 + ^(i)0) I(i)a = (81) 

where 5I^ a is based in g^. 

E. Frame choice 

The energy-momentum tensor in the general (u a ) frame follows from eqs. (21,72,53): 
f ° = - M - J M - ( M + p) V a (V a -B a )-Q a (2V a -B a ), 

T° = (1 - A) [Q a + (ii+p) (V a - B a )\ + ( M + P) (AB a + 2V /3 C a0 ) + + 5p) (V a - B a ) + (V f3 - B fJ ) U a p, 
f$ = (p + 5p) 5$ + 11% + V a [Q{ 3 + ( M + p) (V P - Bp)] + Q a (V - Bp) - 2C or 'IL (h . (82) 

In the energy-frame we set Q a = 0, thus q a = 0. In the normal-frame we have u a = 0, thus from eq. (53) we have 

Energy— frame : Q a = 0, 

Normal - frame : V a - B a + AB a + 2B (3 C a p = 0. (83) 

Although we can take infinitely many different combination of the two frames, the energy- and the normal-frames 
are the ones often used in the literature [22]. By choosing a frame (which is a decision about Q a and V a ) we lose 
no generality. This is because we have 10 independent informations in T a b which can be allocated to the energy 
density fi (one), the pressure^ (one), the anisotropic stress U a p (five, because it is tracefree). The remaining (three) 
informations can be assigned to either the velocity V a (three) or the flux Q a (three); or some combinations of V a and 
Q a with total three informations. 

Thus, in the normal- frame (indicated by a superscript N) we have 
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fg = - fi -6p N -Q Na B a , 
T* = (1-A)Q%, 

if = (p+ Sp N ) 8<$ + TL N % + B a Qp - 2C a m^, (84) 

which is consistent with eq. (78). To the linear-order we notice that S/j,, Sp, Tlap are independent of the frame choice, 
and Q a + (fi + p)(V a — B a ) is a frame-invariant combination [21]. However, to the second-order we no longer have 
such a luxury. As the fluid quantities are defined based on the frame vector as in cq. (23) the values of Sfi, Sp and 
II Q/ 3 are dependent on the frame. 

By comparing the energy-momentum tensor in the normal-frame in eq. (84) with the one in the general frame in 
eq. (82), we find that by replacing the normal- frame fluid quantities to 

Sp N = Sfi+ (V a - B a ) [{fx + p) (V a - B a ) + 2Q a ] , 
S P N = S P +±(V a - B a ) [(p+p) (V a - B a ) + 2Q a ] , 

Qa = Qa + (p + P ) (V a -B a ) + (p + p) (AB a + 2V t3 C a[j ) + (Sp + Sp) (V a - B a ) + (V f3 - B?) U afi , 

= U afj + (n + p) (V a - B a ) (Vp - Bp) + 2Q [a (V fj) - B 0) ) - Igj® {V - W) [(jm + p) (Vy - B 7 ) + 2Q 7 ] , (85) 

we recover the general frame energy-momentum tensor. Thus, by imposing the energy-frame condition (Q a = 0) we 
recover the fluid quantities in the energy-frame. 

Now, similarly, by comparing eq. (82) with the same one evaluated in the energy-frame, we find that by replacing 
the energy-frame fluid quantities (indicated by a superscript E) to 

S[i E = Sfi--±—Q a Q a , 
H + p 

Sp E ^Sp-\-^—Q a Q a , 

6 fl + p 

(p + p) (V E - B a ) = (/i + p) (V a - B a ) + Q a - 2Q?C a p - S -^-^Q a - Q ^, 

t-i + p H + P 



= n a p - -1- [q a Qp ^SQ 7 Q 7 ) , 



(86) 



we recover the general frame energy-momentum tensor. By imposing the normal-frame condition in eq. (83) we 
recover the fluid quantities in the normal-frame. Thus, using eqs. (85,86) we can transform the fluid quantities in one 
frame to the other: 



Sfi N = Sfi E + (v+ P ) (V Ea - B a ) (V E - B a ) , 
Sp" = 5p E + ±(p+ p) (V Ea - B a ) (V E - B a ) , 

Qa = (M + P) (Kf - B a ) + ( M + P) (AB a + 2V E0 C a p) + {5n E + 5p E ) (V E - B a ) + {V E0 - B fi ) U E , 

- nf, + (p + p) (V E Ba) (V E -Bp)- l -g§ (p + p) (V E ^ ET>) (V E - B,) ; (87) 

5fi E = Sfi N l—Q Na QM 

H + p 

5P E = 5P N ~\^—Q Na Q N a, 

3 (j, + p 

(/i + p) (V E B a ) =Qa ~ 2Q N 'C afj - 6 ± + 6P Q% - Q N0^L_ {fl + p) ( ABa + 2B fJ Cap) ^ 

M t P M + P 

nf, = n- - -i- (qZq? - ICQ^Q") • (88) 



F. Spacetime curvatures to the linear-order 

Although straightforward, it is not an easy task to derive the spacetime curvature to the second-order. For our 
purpose, fortunately, it is not necessary to have the forms. Still, it is convenient to have the curvatures to the 
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linear-order and we present them in the following. These follow from eqs. (49,51,52). 
Curvature tensors are: 



aO0 



00/3 



— ) B a , R° aal3 — 0, 



(3) 

9af3 



J\ ' 



-A' + 2 - A 



9 a f3 ~ A ,*\P + B [»\13) + — B Hl3) + C'lp + —C' afj +21—) C a0 , 



a 



R ° a f3~f ~ 2 —9a\J3 A n] ~ B a\[P~<] + TjC^M ~ B (Ml) - 2C 'a 



a -A>5<$ -J*% + \ (Bf + B^)' + \ a - « + B^) + Cf + ±C? , 



-R"o/3 7 — 2— Sf R A^] - B ]R ^ + B a | M ,i - 2 ( — ] 5" R B^] - 2C\ 



\f> A ,i\ - n [p 7 ] +^ |[/3 7 ] 



3 [/3 £. 7] - ZO [(3 | 7] 
2 



w /3 7 a — 



(B 7 I° + B° |7 ) - fl g (B 4 |q + B° |4 ) + 2<5?B (/J | 4) - 2^%, 7) 
g$C? - gfficr + S«C' ps - 5fC' Pl + 2^ - SfC^ 



/3 7 5 



2 a 
a 



(3) 



+^°(/3|5) 7 " ^°(/3|7)5 + °/37 <5 ~ °/3<5 7' 



i?oo — —3 



f >A a -A> + AA- B°( a -B^ a - Cf - a -C«, 

U / U (Ji LL 



i? Q/3 = 2Kg<® + 



- +2 - 



*/3 



(3) R 7 



J./3-" | 7 



+ Ca/3 + 2 ~ < ~'al3 + 2 



+ 2 



Co/3 + ^iX' + 2Cf aWl - - AC Q/3 , 



' / /\ 2 

a 



- ) +if 

a 



6-A' - 12 

a 



- + - 



A-2AA 



+2B a [ a + 6^B« |q + 2CT + 6^-C°' - AKC2 - 2AC" + 2C Q % j. 

The nonvanishing components of the electric and magnetic parts of the Weyl curvature in eq. (20) are: 
E a /3 = — C° 



a0/3 



(l* A \ B \ - \cy - Aq - veer* + <r"^ 



1 , a' 



a 



-,(«)/ 



1 



2 2 



C { $' + (A-2K)C 



7<5 



2 B -y\0)<> + C Ph\s ) - V I 2 ^7l/3)<5 ^Ph\s J > 
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where the symmetrization is only over a and (3 indices. The last steps are evaluated in decomposed forms which will 

9af3 

1 



be introduced in §V. We introduced rj af3 ' r which is based on g^l with rf*®" 1 = r^ a P~<^ and rj 123 = l/\/g^- We have 



~0a/3 7 



a 4 Vl + D 

9 = -a 8 (l + D)g^\ D = 2A + 2C« + 4AC* + B a B a + 2C£C| - 2C£Cf , 



(94) 



which is valid to the second-order in perturbation. 

Deriving the the electric and magnetic parts of Weyl tensor to the second-order using eqs. (19,20) requires quite 
lengthy algebra. Instead, using eqs. (30,33) we can derive them easily. Evaluated in the normal-frame we have 



(I \ (A / CI \ 



B, 17 -B\ + 2C?') 

(a | (a 1 (a J 



K 



ci o 



0h a<*Pft 



\ A — A + -B 1 B^ < 



-a {3) 



^g^g^K aP K j5 + A n (A^ 2C 7 % + C S 5 h ) + A (a A 2 + \bT>b}) 2C s A ms 



H, 



* = { [95% (i - ci) + 2cy 



K, 



(95) 

(96) 



where K and K a p are given in eq. (57); the symmetrization is only over a and (3 indices. We have E a t,u b = = H a i,u b , 
thus E a o = —E a pB^ 1 Eqq = 0, and similarly for H a t. 

For later use, it is convenient to have the spacetime scalar curvature expanded to the second-order. In terms of the 
ADM notation, using cq. (6), we have 



R = R (h) + K al3 K af} + K 2 + - {-K fi + K. a N a - N'- a a ) . 
To the second-order in perturbation, using quantities in §111 B we have 



(97) 



R=R + SR 

(~ 

\a 2 



= 6[^r+H + 2H 2 



cf'V(A + 2/nC Q ° 



-6HA - 12 (H + 2H 2 ) A-2^A + 2 Qb* |q + C°) + 8H Qb* |q + C^j + 2± 
+24HAA - 4A (^B a ]a +C°j - 2 [A + 8HA) Q-B°| Q + C°j + 24 (i? + 2# 2 ) A 2 + ^A^A + ^& a A 
-m l -A^ + ^A Mp C«e + 2^A, (2C^ ]a - C^) + 1 [B {aW B^ + B^ a B^ 2 (B a B°W 
—6HB a B a - 6 (H + 2H 2 ) B a B a 2±B° (2Cf |/3 - C| |a ) ' + 2±B" |a Cg 



An overdot indicates a time derivative with respect to t, with rff = adrj. 



(98) 
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IV. PERTURBED EQUATIONS 



A. Basic equations with general fluids 



In the following we present complete sets of equations valid up to second-order in the perturbation without fixing 
the gauge conditions. As the basic set we consider eqs. (5,8-13,47,48) in the ADM formulation. 
Definition of SK: 



K + 3H 



+SK - 3HA + C° + -B c 



9 • 1 

-A | -HA -C°t- -B l 



\a 



+ -HB a B a 



\ B « (2Cf |/3 - C| |a ) + (c a[i + \B a , 



= An 



where K = K + 5K and K is read from eq. (57) . 
Energy constraint equation: 



167rG/i + 2A - 6H 2 - ^-i? (3) 



+16ttG^ + AH5K 



2Cf '% - 2C Q|/3 



a 13 3 



-SK 2 - 
3 

1 



a 



; "I" -B( a \f3) 



(jap _j_ _J^a\/3 

a 



= Al 

Momentum constraint equation: 



1/3 



- - + -B' 
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+ 2 5K > a + 87rGa( 3c 



= A ( — —SK a - 8nGaQ a ) + A tP 



+ K„-C| l7 )fc + i(s„ h +B\) 



+ 



+ ~ B (a\l] 



+ - 

a 



\a 



= N 2a . 

Trace of the ADM propagation equation: 
- 3H + 3H 2 + AttG {(i + 3p) - A 

+SK + 2HSK - AttG {Sfi + 35p) + ^A^ a a + 3HA 
= A5K - \5K, a B a + hx 2 + ^H (3A 2 - B a B a ) 



(99) 



(100) 



(101) 
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+— 

a z 



2AA\ a a + A a A' a - B?B f } a a - B^ a B fj[a + A' a (2Cf |/3 - C| |a ) + 2C Q/3 A 



"1/3 



= JV 3 . 

Tracefree ADM propagation equation: 



(102) 



3^ 
1 



|a 



3iJ 



C 



a 



17 



C7 + -B J , +3H C1 + -B 



1 ' 2,/ 
1 



/3 



?7 



a 



17 



C" 



+ ci?'\ cf\ c:} a R Ir^cs - h a a ( 2ct h K - 2c:}\- 2 -R^a 



|/?7 T ^i 9 7 W /3 7 ^7/3 3 Jt ' ^P 3 U /5 1^7 5 ^7 <5 3' 



87rGn^ 



2a 



A + 2C" 7 ( C f)l "I" -£>, ; , 



B (/ 3| 7) j +-^ 7 (c a7 |/3 + C ( 



7|a ^~<a\~f 
& °/3 



+ 3i? 



r<a7 1 r>"y\ a r< a h 



1/3 



a 



<S? + 5(*-„ + V) 



B~< + 5K 



|7 



^ + 2^ ( B I 



1/3 ^ S /3 



1 



+3F 



2C-A/3|7-(^ 7 |, + ^ 7|Q -^ 17 ) A7 



CI. + --B 7 i A + 2C 75 C 75 + -B. 



B d I 2C] h - C 7 5 



+ [CI + -B 7 



A 



1 



C-< + -B^ b \ B d +5K[C: + -B 



1 



|7 



+ - 



a I/ 3 



^ + l(fl- |7 + B 7 l-)_ 



-Is?, 



ir • f q |5e + c 7|e5 - c 7 , £ e - q 75 ) - -r^c*cj 



r (2^,, - q 5 ) (2c 



-,7(5 



|7 7 



7<5|e 



-167rGC a7 n^ 



4/3- 

Energy conservation equation: 
fi + 3H (fi+p) 



(103) 



+5(i + 3H {5fi + 5 P ) -(ji+p) (5K - SHA) + -Q a ]a 
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= --S^ a B a + {Sfi + Sp) (SK - 3HA) + (pi+p) 



ASK + -H (A 2 - B a B a ) 



AQ a {a + Q a (2A. a + C| |a - 2Cf |/3 ) - 2C afi Q aW 



n Q/3 ( C a(] + -B a \ P 



= N 5 . 

Momentum conservation equation: 
1 
a 



(104) 



Q a + 4HQ a + ^- {fjL+p)A ia + 5p, a +IL^ 



Qa (SK - 3HA) + - 1 - Q*\pB fj - QpB P w - (5fi + Sp) A, a + A [( M + p) A a - Sp, 

- ( M + P ) B l3 Bj } \ a + 2 (c^n Q7 ) |/3 - n^ |7 + n^c^ - 



n 



a|/3 



6 a ■ 



(105) 



In the multi-component situation we additionally have the energy and the momentum conservation of individual 
component in eqs. (47,48). 

Energy conservation equation for the i-th component: 
A(») + 3-ff (//(i) + + --f(»)o 

+ 3iJ + - (m( 4 ) (SK - 3HA) + ^Qf i)la + \ 5I (i)o 

= -U^ a B a + (S m + 5 P{1) ) (SK - 3HA) + ( m + p {i) ) ASK + ( m + p (i) ) (A 2 - B a B a ) 



1 

a 

(i)5- 



^[- a B a \ [} + C a p\-- a 5I {l)a B c 



1 



1 



(106) 



Momentum conservation equation for the i-th component: 

<9(»)q + ^HQ (i)a + i + p {€) ) A, a + i ^P(i) )Q + n ( . } f |/3 - W( i)Q ,) 

= a | (^W'" + n (»)"l/3 ~~ SI (*)<*) ^ ~ (<^« + 5 P(i)) A« + (M(i) + (^A<* - B f) B p \ a ) 

-Q{i)a\(iB (i - Q {i)P B^ a + a(SK- 3HA) Q (i)a + 2 (C^n (i)£ry ) |/3 - C|, 7 II ( + C^II^ - A, p U ( / a j 

= N {i)6a . (107) 

The collective fluid quantities are given in eq. (75). The equations are presented with the quadratic combination of 
the linear order terms located in the right-hand-side. Still, notice that the equations are presented to the second-order 
without separating the background order part. 

Equations (99-107) provide a complete set valid for the Einstein's gravity with an imperfect fluid, thus the most 
general form of energy-momentum tensor. We have not imposed any condition like the gauge condition. In the 
following subsections we will consider the cases of minimally coupled scalar fields, an electromagnetic field, and a 
broad class of generalized gravity theories. We emphasize that even in these additional fields or generalized gravity 
the above equations remain valid with the fluid quantities reinterpreted to absorb the contributions from the fields 
and the generalized gravity. 



B. Scalar field 

1. Covariant equations 
The action for a minimally coupled scalar field is given as 
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d 4 x. 



(108) 



The equation of motion follows from the variation in (f> 

& C c-Vi = 0, 



(109) 



where V; 



^j. From Sg ab CM = hV—gT ab 6g a b we have the energy- momentum tensor 



f { a i ] = 4>J, b - -~g ab 4>' c h - ~g ab V{4>). 



(110) 



2. Perturbations 



We decompose 

4> = <f> + 5(f). 

The equation of motion becomes 
4> + 3H<j) + V4, 

+S4> + 3H6<j> - ^5<f> + V^Scj) -2A4>-<j)(A + 6HA - -B a , - C% 



2A 



(111) 



+<f)B a B a + 
= iV . 

The energy-momentum tensor gives: 



6<j> + 3HS<j) -2A4>-j> \2A + 6HA - -B a ]a - C°) +6<j)[A- -B a ]a - C°) - 2—B a 8(f) . 

v (x J \ Ob /a 

-5<f>, a (-A< a — B a — 2HB a - 2-C a ? + -C^ a ) - 2^80 aW C a0 - ]-V 4,4,48 ft 
a \a a |p a p / a z ' 2 

~^A a B a + B a B a + 3HB a B a + ±B fj (2C Q/3 |a - C^) + 2C af) (^B a \ fj + C af: 



(112) 



1 oo 



a 2 V + ft 5ft + a 2 {V ^(f, + 2VA) 



+ \a 2 V 448ft + 2a 2 V^5<pA - ft8(f,. a B a + X -ft 2 B a B a , 



T, 



(4>) _ 



TW _ „(3) 
1 afJ 



'H,a - {^ft 2 - a 2 V^j B a + 6<f>'5<t>, a + {-ft 8ft + a 2 V^8(f) + ft 2 A) B a , 

9 ( £ (^' 2 - « 2 ^) + 9% ( W - a 2 V^54> - ft 2 A) + (ft 2 2a 2 V) C a(i 
+54>, a 5<t>,p + 2 [ft 5ft - a 2 V^ - ft 2 A] C af3 



~\9«l M^' 2 + W^W 1 + a 2 V^5ft + Aft 5ft A - 2ft 5<f)^ + ft 2 (-AA 2 + B 7 B 7 )] . 
Fluid quantities can be read from eq. (79) as: 

+ SftW = ^ft + V + <}>5<f> - ft A + V 4.5(f) 

+ \sft + ^S(f), a S(f)' a + \v 4,4,5ft - 2<j)5j>A + ^ft<P, a B a + 2ft 'A 2 - l -ftB a B a 
p W + s p W = ift _ v + <f>S<f> - ft A 
+\sft - ^6(j) ta 6(j)' a 



(113) 



\v 4,4,8ft - 248(f) A + -484>, a B a + 2ft A 2 - \ftB a B a 

A Qj A 
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+ 



5<p, a (5<j) - (f)Aj 



,(3); 



(114) 



We indicate the quadratic parts (the quadratic combinations of two linear-order terms) as 5jjS q \ 5p^ q \ Qs , and Il^g. 



C. Scalar fields 

i. Covariant equations 
The action for multiple components of minimally coupled scalar fields is 



-<) 



1 ^ 1 



where = 1, 2 ... n indicate the n scalar fields. The equation of motion for i-th. component is 



vr -vj =0. 

(») c ,0(i) 



The energy-momentum tensor is 



^if = XI ( ^{k),a<l>(k),b - ^9ab<t>( k f<l>{k),c ) ~ <?a6^(</>(0 )• 



(115) 



(116) 



(117) 



2. Perturbations 

We introduce 

= 0(i) + 

The equation of motion for i-th component becomes 



2A 



i)\A--B a \ a -C%\-2-B a 5<t> (i)ia 



(118) 



+ -#(*),<* 



+0 (l) B Q B a + (4) 
The energy-momentum tensor gives: 



' k,l 

\A, a B a + B a B a + 3HB a B a + I Bp (2C af \ a - C^) + 2C a ? (^B a]p + C a[j 



(119) 



k k k \ 

+ \ a2 Yl ^m^o^W^W + 2a2 ^ (fc )^(fe) A - <t>'(k) S< t>(k), a B a + )-4>f k) B a B a J , 
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(k) 0< P(k),c 



\ E - B - + E [^(fc)^(*),a + (-4>[ 

\ k Ik 



= flS E - a2 M + ffS E 0'w<%) - - €) A ) + ( E - 2a2y ) c ^ 



+E 

+« 2 E^ 



(k),a°<P(k),0 



+ 2 



fc) 



*(fc) + ^(fc),7 5 ^(fc)' 7 



(fc) <ty (I) + 4^ (fc) <ty' (fc) i4 - 2^ (fe) <^ (fe);7 iF + 0f fe) (-4A 2 + £ 7 iF) 



(120) 



Fluid quantities can be read from eq. (79) as: 



M W + = ± £ 4?^ +V + J2 (W^(fc) - % k) A + V,+ W 5<j> w ) + E 



k L 



+i E ^w^o^w^o - 2 ^(k)S>i>(k)A + -4>(k)S(f>(k), a B a + (2A 2 - \-B a B a \ 4> 2 (k) 



P w + 5 P w - e 2^) - v + E (^(fc)^(fc) - - *W)<^w) + E 

k k k 

"2 E^«= 

= -~E [^(fc)^(fc),a + (<*0(k) -^(fc)- 4 ) 50(fc), 



^(fc) - 6^2^(fc),a^(fc)' 

2 ^ ■ , vl H)<i>m S< t>(k)8<t>(.i) - 24>(k)H(k)A + ^<j)(k)$4>(k), a B a + (2A 2 - l;B a B a j <^ fe) 



We indicate the quadratic parts as <5^ 9 ) , , , and 11^ 



a/3- 



(121) 



D. Generalized gravity theories 



1. Covariant equations 



As the action for a class of generalized gravity theories we consider 



V(<j> K ) + L r , 



d 4 a 



(122) 



(j) 1 is the 7-th component of TV scalar fields. The capital indices J, J, K, . . . = 1, 2, 3, . . . , N indicate the scalar fields, 
and the summation convention is used for repeated indices. f((j> K ,R) is a general algebraic function of R and the 
scalar fields (j) 1 , and gij((j) K ) and V(<j) K ) are general algebraic functions of the scalar fields. We include a nonlinear 
sigma type kinetic term where the kinetic matrix gjj is considered as a Riemannian metric on the manifold with the 
coordinates c/r . The matter part Lagrangian L m includes the fluids, the kinetic components, and the interaction with 
the fields, as well. We have introduced the general action in cq. (122) in [23,12] as a toy model which allows quite 
general handling of various different generalized gravity theories in a unified manner, see [24]. Our generalized gravity 
includes as subset: f(R) gravity which includes R 2 gravity, the scalar-tensor theory which includes the Jordan-Brans- 
Dickc theory, the non-minimally coupled scalar field, the induced gravity, the low-energy effective action of string 
theory, etc, and various combinations of such gravity theories with additional multiple fields and fluids. It does not, 
however, include higher-derivative theories with terms like R ab R a i,, see [25] for its role. 
The gravitational field equation and the equation of motion become: 
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G ab = j 



T„,, f <!,.,[ ft, a ft, b - \g a bft' c ftA +l(f-RF- 21 ) + i.. rJ , - ~g ab F" 



(123) 

(124) 
(125) 



where F = df/dR; g IJ is the inverse metric of g u , = \g IL (§lj,k + 9lk,j - 9jk,l), and Vj = dV/(dft). 

Introduction of the effective energy-momentum tensor Tffjf^ provides a useful trick to derive and handle the perturbed 
set of equations [23]. It allows the equations derived in Einstein gravity remain valid with the energy- momentum 
parts replaced by the effective ones. 

We note that the gravity theory in eq. (122) can be transformed to Einstein's gravity through a conformal rescaling 
of the metric and rescaling of one of the fields. As the result we have Einstein's gravity sector with only complications 
appearing in the modified form of the field potential; the nonlinear sigma type couplings in the kinetic part also 
remain. We have studied the conformal transformation properties to the linear order perturbation in [23,26], and in 
most general form in the Appendix A of [29]. Extention to the second-order perturbation is trivial. 



2. Perturbed equations 

The perturbed set of equations can be derived similarly as in the previous sections on the scalar fields. We set 



Thus, 



f = f + 5f, f J = r J + £r 7 . 



SF = Fjdft + F R SR + ^FjjSftdft + F JR 5ftSR + ^F >RR SR 2 . 



(126) 



(127) 



= 2A 



'. f~iOL 

\a U » 



2-B c 

a 



2—C a 



,a|/3 



The equation of motion in cq. (124) gives 

ft + 3H~ft - l -g IJ (f,j - 2Vj) + Y I JK ftft i + T 1 

+5ft + 3HSft - ^5 ft - ^Aft + ft (-A - 6HA + -B a {a + C°) - \g IJ [Fj5R + (f, LJ - 2V LJ ) 6 ft] 
if, j - 2Vj) + 2Y I JK (ft 5ft 1 - A<j> J <f> K ) + T^Sftft^ + 8T 1 
6ft + 3H6ft - 2Aft - ft (2A + 6HA - ^B a {a - C^j + 2T I JK (j> J 5<j> K - A ft ft 1 *} + Tj KL 8ftft(j) K 

+B a B a (ft + wft + r^ft^} + (a-- 

a ( ~A' a — B a — 2HB a - 2-C af \ R + -Cf 1 
a ' \a a 1 ' a ' 

\A, a B a + B a B a + - a B a (2C a ^ - C^ a ) + 2C^ (^-B a \ (j + C afj 

[F.rjSR 2 + 2F LJ 6R6ft + (f, LM j - 2V LM j) 8ft 6^ M ] 
+ \d 1J jft [FjSR + (l L j - 2V LJ ) 5ft] + l g IJ LM 5ft50 M (f,j - 2Vj) 

+Vjk (-Sft6j> K - 2-ft8$ K a B a + Isft^d^ - 2Y l JK ^ftftbft l - V 7i ^ 
= N g . 

SR can be read from eq. (98). From eq. (123) the effective energy momentum tensor gives: 



+4>- 
1 



(128) 
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8irGF 



Too + \gij¥'$" - \a 2 (/ -RF- TV) - 3^F' 



+g IJ <P I '54> J ' + \gu^4> L ¥'4> J ' - \* 2 [(f,L - 2V L ) 54> L - RSF] 
-a 2 A (f -RF- 2V) - 3—SF' + ASF - (b\ + C£") F' 

+9IJ 



- A gu,LMH L ^ M ^'^' 



+ Ja 2 [F^SR 2 - {/,lm - 2V LM ) 6^ L S(f> M ] - a 2 A [(/, L - 2V L ) 6<b L - R6F] 

- (B a \ a + Cf) SF' - 2B a SF' a - (^-B a + 2C Q/3 |(3 - SF. a + 2AASF - 2C af3 SF MP 

2A, a B<* + B a (2C a % - Cff a ) + 2<7°" J (B a p + C aP )] 



+B a B a ( F" + —F' ] + 
a 



i 



8nGF 



f aa + gu^'H^a - 7;B a [gutfW + a 2 (/ - RF - 2V)] + 6F[ a - -5F >C 



-B a F" - \A, a + -B a j F' 

+ 9IJ [Stf'S^^ + B a (-^'Scf,-" + W)] + 9lJ ^ L (■ 



2 Ba ' 



\a 2 B a - 2V L ) H L - R6F] - B a SF" - (a i<x + ^B a ) SF' - (\b^ - l -B ( \ a + C%) SF 



+B a ASF + 2AB a F" + 



2AA, a + 2^-AB a + BpC* - B f3 B [fj{a] + B a (a 1 - B^ - C% 



F' 



gu<t> 4* + a 2 if -RF- 2V) + 2F" + 2-F' 



+o (3) 



gu4> 4* + a 2 (/ -RF- 2V) + 2F" + 2-F 1 



+ SF ,a\P ~ ( B (a\l3) + C' a p) F ' 



g u {^'S^' - Atf'tf') + ^u,Lty W + \a 2 [(f,L - Wl) Scf> L - R6F] 



+SF" + -SF' - ASF - 2AF" — [ A' + 2— A - B 1 , -C!'\F 



+ gi.jS4> 1 • ' 



+c, 



a/3 



2gu {4> D 54> JI - A^'cj> J ') + gurf^tfW + a 2 [(/, L - 2V L ) S<j> L - RSF 



+2 



SF" + —SF' - ASF - 2AF" - [ A' + 2- A - B 1 ^ - C% ) F' 



(B (aW + C' ap ) (SF' - 2AF') (2Cj am - cj") (SF n B, ( F') 



+a (3) 



-gu [Stf'Scj,- 1 ' - AA^'S^' + {AA 2 - B^B^) + 2B^ V S<\>\ 7 - Stf^Stf ' ] 



+9i.j,lH L (0 



W J ') + - A gu,LMH L 



+ -a 2 [-F iR SR 2 + (l LM - 2V LM ) 8<f> 
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-2ASF" - ( A' + 2- A - B~< h - CI] ) SF' + 2B^5F' 1 



+ \-A« + B 1 ' + + 2C' ,s ls - C d s h j 5F„ + 2(P S 6F„\ S + (4A 2 - B^'B y ) \F" + --F' 
AAA 1 - A a B"< - STB; - 2A (V, 7 + G>') - S 7 (2C^ ls - C 5 5 ' 7 ) - 2C~< 5 (B y{s + C' yl 



+ 

The fluid quantities follow from eq. (79): 
1 



F' 



(129) 



/i (eff)+^(eff) 



^ F ^+\guW-\{f-RF-2V)-*HF 

+ g u (pStf - A<j> T <j> J ) + ^g IJlL S4> L ^4> J - \ (fx - W L )5ct> L 
-3H5F + (\r +^)SF- (-6HA + V, Q + C°) F 



:9u 



^54> J + -l<ty I,a ty J , Q - AAtfSp + 2-B a <j> I 6cf> J ta + (4A 2 - B a B a ) ft j> J 



+9i.j,l5<P l (<W j - A4> T 4> J ^ + ^guxuH 



L XJ.M1I1J 



+ - [F R 5R 2 - (f, LM - 2V LM ) 5<t> 



3HB a + i (2C afi ](j - C^ |Q ) 



- 6HA+ -B\ Y + C* 



SF >a - B, 
a 



(5F - 2AF^j 

■ a F) - 2^C^SF M0 + 2FC a < 1 (^B aW + C af ^J j , 



,w + Sp w = 



■^PY + \9™V* 3 + \{f~RF-2V)+F + 2HF 

+S P + gu (<W J - + \gu, L H L ^^ + \ (f.L - 2V L ) 5^ L 



+SF + 2H5F - ( -R + ~ ] 5F - 2AF 

2 6CL Z 



A + AHA- 2 -(^B\ t + CZ 



+9u 



^5<t> I 5<}> J - 2A^5^ J + i (4A 2 - B a B a ) ft 4> 3 + ^B a ft 5<p J , 



1 

6^' 



+9i.j,lH L (ftsft - Aftft) + X -gu,LM^ L ^ M ¥tf - \ [F, R 6R 2 - (f, LM - 2V LM ) H 

A + AHA -%( -B a , + C°) (SF - 2AF] + 2-B a 5F a 
3 \a J J V / a 

\A<" + B"+HB° + ±. (2C a( \ p Cf a )j (-JF, a - B a F^j + ^C af) 5F Mf3 



-2 ASF 



+ (4A 2 - B a B a ) F + 



2AA - -A a B a - HB a B a 

a 3 



C af3 { -B a \p + C al 3 



F 



SirGF a 



1 



a 



-gu 



5F a +H5F c 



I 1 - A^ 5<f> J ia - - gi j, L 54> L ¥5ct>\ a + A 



-A a F 
a 

-3-FA a + - (sF a - HSF a ) 



+ -A. a SF+- 

a a 



n (eff) = 



Q/3 8irGF \ 



I (3) 



8F, a \p ~ [ - B (a\p) + Co-Si ) F - -g t 



*SF- (-if . ■ c: F 
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1 T j 2 



;5F 



B\ + C])F 



(±B (aW + C a ^j (6F 2AF) \ (2C? aW cj~<) (±<JF i7 - B,F 



(3) 



\giM^84> J a + 2C~< S (-B 7]s + C jS ] F - ^ S SF^ S 



(130) 



We have used eq. (78) for the energy-momentum tensor. We indicate the quadratic parts as 5fi( eS ' q \ Sp( cS > q \ Qa S ' 9 ^ 



and n^' 9 '. We note again that in this generalized gravity the basic equations in §IVA remain valid with the fluid 
quantities replaced by the effective ones. 



E. Electromagnetic field 

1. Covariant equations 



The Lagrangian of electromagnetic field is given as 



(131) 



where F ab = A a , b — A b ^ a . The energy-momentum tensor becomes 



p(em) 
L ab 



1 



F ac F b c - -^9abFcdF C . 



We introduce [20] 



Then we have 



^(em) 



-(em) 



F ab = u a E b - u b E a - r) abcd u c H d , 

E a = F ab u\ H a =^f) abcd u b F cd , E 2 = E a E a , H 2 = H a H a , 
q = -j a u a , J a = hfj b . 



^u a u b (e 2 + ff 2 ) + 2u (a f, b)cgd u c E 9 H d - E a E b - H a H b + H ab [e 2 + H 2 ) , 



\[e 2 + h 2 



) 



p(em) 



= \{e 2 + h 2 



)■ 



(em) 



lr 



-E a E b — H a H b + g^ab ( E 2 + H 2 



From the Maxwell equations and the conservation equations 

F ab ., b = ~f, F [ab , c] =Q, ~f. a = 0, 
we can derive the covariant forms of relativistic Maxwell's equations [20] : 

E a . b h b + 2H a u a = q, 

H a . b h b a - 2E a Cu a = 0, 



h a b E b :c u c = E b (q\ + ~a\ - l§h$) + f) abcd u b (a c H d - H c , d ) - J a , 



KH b ., c u c 



H b (u, a b + ~o\ - lfh a b \ + fj abcd u b (a c E d - E c , d ) - J a , 



q, a u a + eq + h a b J b . a + J a d a =0. 



(132) 



(133) 



(134) 
(135) 

(136) 
(137) 

(138) 

(139) 
(140) 
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2. Perturbations 



We take the normal-frame, thus u a = n a , thus uj a b = 0. Due to the high symmetry the Friedmann background does 
not support the electric or magnetic field. Thus, E a and H a are already perturbed order. We set 



TP — TP ff = ff 



(141) 



where E a and H a are based on gf^. Thus, E = -E a B a (which follows from E a h a = 0), etc. For rf bcd see eq. (94). 
Equations (136-140) become: 



E« la -aH q = 2(C^E fj ) la -C: w EP, 
H a la = 2(C^H ) la -C^H", 

E a + HE a + \ aP ~>H pb + J a = A^E a + HE a } - -\e\ } B p + E Q-B"^ + 2C^j 



5q + 3H6q = -3HASq - -5q a B a + 5KSq - \ 

a a-* 



{r-2c^j fj ) ]a + r (cl ]a + A, a ) 



where we set 



q = q + Sq, J a = J a , 



(142) 
(143) 



(144) 
(145) 
(146) 

(147) 



with J a based on gf^; J a in this subsection differs from the flux term in ADM notation used in the other sections. 
We have q = 0. 

The energy-momentum tensor becomes: 



n(em) 



1 0a --VafhEPir, 

TiT = -E a Ep H a Hp + IgW {E^E, + H~<H,) . 



Fluid quantities can be read from eq. (79) as: 



5^ = 35^1 = JL (E a E a + H a H a ) , 



n 



(cm) 
a/3 



E.E « + /./„ U « - p<® (WE y + HiH n ) 



We have M (om) = = p<- cm K 



F. Null-geodesic and temperature anisotropy 



(148) 



(149) 



We introduce the photon four-velocity as: 

k° = -(v + 6v), k° 
a 



{e a + 5e a ) : 



1 + — + 2A- B a e a + 2A— - B a 5e c 
v v 

k a = -av (e a + Se a + B a + 2C al3 e f3 + B a —j- + 2C a/3 Se^ , 



(150) 
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where e a and Se a are based on g^p, and v and e a are assumed to be the background order. We have 



d dx a d 



Sv. 



Thus, 



-jr = ~~^T~^~~^ — ka da = -[do- e a d a + —8 - Se a d, 
dX o\ ox r 



— = d - e a d a , 
ay 



(151) 



(152) 



can be considered as a derivative along the background photon four-velocity. The null equation, k a k a = 0, gives 



k a k a = v 2 



Sv 



e a e a - 1 + 2 ( e a 5e a -^--A + B a e a + C a0 e a e 



X 2 x 

+Se a Se a -JL-2—(2A- B a e a ) + 2(B a + 2C« /3 e /3 ) Se a 



0. 



The geodesic equation, k a . b k b = 0, using eq. (52), gives 
v 2 J (av)' ( Sv\ v' Sv Sv 



k^ .fok 



a* av 



+ 2- 



v v 



e a + 2-e a Se a + A' - 2- A 
v a a 



+ (B a]0 + C' ap + ^C a ^j e a e fj - 2 (a^ - ^-B^j e a 



Sv Sv' Sv 



8v Sv 



'-Se a + 2— A' - 2— A, a B a e a + —Se a 6e a 



~25e a (A ia - ^b}j - ^e a Se a A + 2e a 5e> 3 (c' af3 + 2^C afj + B H/J) ) 

A (2B alP + 2C' afJ + ^C af3 \ + B, (2Cl w - cj 7 )] + £ (4A 2 - B a B a ) 



~k a A b = -e°" + e e a ,» - Se a ' - —e al + Se a u ,e p + 8e t3 e a , 



1/3 



\0 l 



1/3 



+ (2C£| 7 - cj") tPe* - (B f } a - B\ + 2Cf) e (i + A- a - B al 

~Se a ' + 2^ (A a - B al ) - (de? + (b^ - B« |/3 + 2Cf) + 8eHe<* |j9 

+2e t3 5e< (2C£, 7 - cj") + A'B a - 2A. fj C afi + 2C%B fj ' - 2B a A^ + 4C a ^ (B [m + C 

-2C<i (2C% y - C^ S ) eV + B a (B 0h + C' M ) eV j = 0, 



where we used the null equation in eq. (153). To the background order eqs. (153-155) give 

V 



e a e a = l, v(xa-\ e a ' = eV. 



Using eqs. (152,153), eq. (154) becomes 
d f Sv 



dy \v 



+ A\- A. a e a + (B alfJ + C' a0 ) e a e' 3 



Sv Sv' a' Sv 2 Sv 



v v a v* 



Sv . 



Sv 



-2—\A' + 2— A + ^^Se a + 2— A 



(153) 



(154) 



(155) 



(156) 
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+2Se a A. a + 4-e a 6e a A - 2e a Se' 3 (B (aW) + C' afj ) 



2A B, 



C'a/3 + ^Cafj 



e <* e P _ ± ( 4j 4 2 - B a B a ) 



-2 (^2AA a - 2^-AB a + B p C p a ' + B fj B [am ^j e a + 2AA' + A. a B a - B a + ^-B c 



N v . 



(157) 



Thus, we have 



8v 
— + A 
v 



[° [A. a e a - {B alfj + C' af3 ) e a e? + N v ] dy, 

J E 



(158) 



where the integral is along the ray's null-geodesic path from E the emitted event at the intersection of the ray and 
the last scattering surface to O the observed event here and now. 

The temperatures of the CMB at two different points (O and E) along a single null-geodesic ray in a given 
observational direction are [27,17], 



To _ 1 (~k a u a )c 



f E l+~z {k b u b ) E ' 



(159) 



where u a at O and E are the local four-velocities of the observer and the emitter, respectively. Thus, u a should be 
considered as the one based on the energy-frame which sets q a = 0; or equivalently in a general frame vector which 
absorbs the flux term to the frame vector to the second-order. 
Using cqs. (53,150) we have 



k a u a = -v 



l + -+A+(V a E - B a ) e" + -A + 6e« (V a E B a ) + (AB a + 2C a(j V^) e<* + \v Ea V E \tf 



(160) 



We have denoted the energy-frame nature by replacing V a to V E ; if we consider eq. (86), eq. (160) in this form is 
valid in the general frame. Since the calculations in the rest of this paper are based on the normal-frame vector we 
use eq. (88) to derive the result in the normal-frame. We have 



k a u a = -v 



1 + ^ + A+ J_ Q a _ (S, + 5 P )Q a +U aP Q^ 



+ ^A + 5e«- Qa 



H + p 



Qa 

fi + p 



B a + 



Q° 



n + p 



2 



(161) 



Using eq. (159) we have 
T E vo 
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2 



Sv 

r A -\ 

v H + P 
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Q- ! • y 



A 



H + P 



N 7 



(162) 



where = 4j£ 



= ^.\ - *f-\ E and %f\ E = %f at E. Thus, if we take T Q /T E = v Q /v E , eqs. (162,158) give 
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ST 



ST 
~T 



1 



r° 

/ [A, a e a - (B alp + C' af} ) e a eP + N v ] dy + N 7 

J E 



(163) 



In the large angular scale we are considering (larger than the horizon size at the last scattering era) the detailed 
dynamics at last scattering is not important. The physical processes of last scattering are important in the small 
angular scale where we need to solve the Boltzmann equation for the photon distribution function, see §IVG. 



G. Boltzmann equation 

1. Covariant equations 



The relativistic Boltzmann equation is [28] 



d_~ d_x^df_ dj^dl df_ _- h df_ - ~ 

d\ J dx dx a d\ dp a p dx a bcPP dp a uu 



(164) 



where f(x a ,p b ) is a distribution function with the phase space variables x a and p a = dx a /dX, and C[f] is the collision 
term. The energy-momentum tensor of the collisionless (or collisional) component is 



f£>= [ 2e(p°)S(p% I nr u>., Pl J- ' 



.9 



28(f)S(p c p c + m')p a p b f^-~gd 4 p^, 



where 



r<~r~ 9, r/~ -n ~n n ~ -n -r, ~ -r. ~r ?n <5(mass shell) 5(mass shell) 
S{p c Pc + to 2 ) = S(g 00 p p + 2g 0a p°p a + g a fiP a p^ + ™ 2 ) = ^ '— = -i '- 



\2~g m p* + 2go a p° 



2|p | 



Thus, after integrating over p°, we have 



-gd 3 p 



-123 



-PaPbf, 



(165) 



(166) 



(167) 



with the mass-shell condition p a p a + m = 0. 

2. Perturbed equations 
Under our metric, using p a as the phase space variable, we have 

P°f + P a f, a | ^ (PV + 9 { $P a f) + A'pV + 2 (A, a a -B}j fp« 

+ (-4 3 S a + b «\p + c u + 4^) rf 

+ (-2AA' - A, a B a + B a B' a + ^B a B a ^j p°p° + 2 (-2AA, a + 2^-AB a - BpC?' + B^B Wc ^j p°p a 



+ 



p a f 



(2^fp a + T (3) ; 7 /p^ + (a* B a > ^B a ^j p°p° + (B^ a - B a {0 + 2Cf) pY 



df 
dp 



S-gf^ + 2C^-C^)^ 



df_ 

dp a 



C[fl 



(168) 



As the phase space variable it is convenient to use (q, 7") introduced as 
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q = a\J a 2 (jp) 2 — m 2 



yVOp ) 2 - to 2 

a 2 (p ) 2 
a 2 (p°) 2 — m 2 



a 2 (p ) 2 



1 + 



p a 



a 2 (p ) 2 — m 2 
a 2 {p ) 2 

n 2 (ffi\ 2 — m 2 



A+-B a B a 



1 a 4 (p ) 4 



2 [a 2 {p ) 2 -m 2 } 2 



A 



Bp +C,p + 2[a2{p0)2 _ m2]2 A P 



+ 



AB a f - AC^f - '-B Si B (i p a ) + C%B f3 p° - ^C}C«f 



P° = \^q 2 +m 2 a 2 (l - A + %A 2 - l -B a B a 



a" • ~ \ 2 2 

\ qi a + vV + m 2 a 2 B a - qC^ 3 - V ' q 2 + m 2 a 2 (AB a + 2C%B ) + \qC}C^ 

where 7 a is based on g^. The mass-shell condition gives 

(P ) 2 - g%P a f - m 2 /a 2 + 2A(p ) 2 + 2B a p a f - 2C a pp a f = 0, 
and we can show j a j a — 1. The Boltzmann equation becomes 



/' + /2 q 22 ( rsf, a - r (3) ; 7 7Vf^) 
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a 2 J q 



\f q 2 + m 2 a 2 



+ 



\J q 2 + m 2 a 2 



l + A-±A 2 + ±B a B a )(lf} 



For convenience we located the collision term in N c . The energy- momentum tensor becomes 



1 ab 



Pbf 



q 2 dqdfl q 
sjq 2 + m 2 a 2 



Thus, using / 7a jpdQ q = §3^, we have 
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From eq. (79) the fluid quantities become: 

M ( c > + 5^ = ^J vV + m^q 2 dqdn q (/ + Sf) , 



3a 4 J ^qi + m 2 a 2 

Qi c) = -A q 3 dqdn q ( 7a + c a/J y») 



If we have multiple components each described by the Boltzmann equation, all equations in this subsection remain 
valid for any component with / replaced by f^, etc., and the total (collective) fluid quantities as the sum of the 
individual ones. 



V. DECOMPOSITION 
A. Three perturbation types 



We decompose the perturbation variables as follow: 

A = a, 
B a = 0, a + Bg\ 



n a/3 = 1 (n \ g % An) + inj^ + ng, (175) 



with the properties: 



B [v)a , =0, c (v)a , =0, v {v)a , =0, n (u)Q , =0, 

I a ' I a 1 I a ' \a ' 



Cr(t)a= 0) n W a a = 0, 0^ = 0, n w * a|/J = 0. (176) 

/i + p appearing in the decomposition of Qq, is assumed to be the background order quantity. The vector- and the 
tensor-type perturbations are denoted by superscripts (v) and (t), respectively. We assume all these variables are 
based on g^. The decomposed variables also can be expressed in terms of the original variables. For example, we 
have (3 = A _1 V Q i3 Q and Ba ^ = B a — V a A _1 V^.B/3, etc, where V a means vl 3 '. For the fluid quantities we have 



q { : ] = Qa~ v q a- i v 3 o , 



II |r (A • ) A \ \ II 



ng = *<# - \ (v-vp - \ g § a) (a + _1 a- 1 v Vn* 

-2V (Q + |i? (3) ) (V^n % - V^A-^V^f) • (177) 



We introduce 
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(178) 



and let 

k = SK. (179) 
In the multi-component situation we have eq. (75). For the individual component we have 

Q(i)a = (P(i) +P(i)) (-"(»),« + ' 

,(3) att. ^ _l IrrW xnW 



n (i)a/3 ^ ^ ^n ( o,ci/3 - 3^An (i) j + -n$ (a|/J) + n$ a/J) (iso) 

with 

' 7 (i) |a = u ' |a = u ' 11 (i)a = u ' 11 (i)a|/3 



^ )Q , =o, n<" )a , =o, n^ )a = o, n^ 3 ,,, = o, (isi) 

fx) a ' m a ' ma ' u a d ' V / 



and 



<ji (0a = *j (0 , a + jj$ Q) «5/£ )Q |Q = o. (182) 



The definitions for the scalar-type perturbation variables are introduced to match with our notation used in the linear 
analysis [7,12,29]; compared with our previous definitions in the linear theory our v and correspond to v/k and 
i>(i)/fc in [29] where k is the wave number. These are the notations introduced by Bardeen in 1988 [7]. A complete 
set of equations written separately for the three perturbation types will be presented in equations (195)-(210) where 
the quadratic combinations of the linear-order variables contribute to the second-order perturbations. Thus, to the 
second-order the three perturbation types couple each other through quadratic combinations of the linear-order terms. 
If needed we may decompose the perturbed order quantities explicitly as in eq. (50) 

a = a^+a^\ <p = ^ + (183) 

etc. 



B. Background equations 



To the background order, eqs. (100,102,104,112) give 



, 8nG K A , 

H =—»-^ + r (184) 

a 4nG „ , A . . 

-=-—(M + 3p)+3. (185) 

/i + 3H(M + p) = 0, (186) 

$ + 3H4> + V,4, = 0. (187) 

In the multi-component situations from eqs. (106,119) we have 

A(») + (m« +?>(»)) = ~-I{i)o, (188) 

+ 3Hj> {i) + V t4>w = 0, (189) 

with 

k k 

which follow from eq. (121). 

In the generalized gravity considered in §IVD, eqs. (184,185) remain valid by replacing the fluid quantities to the 
effective one in eq. (130): 
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(off) 



8ttGF 
1 



(/ -RF- 2V) - 3HF 



P + \gijW + \ (/ - RF - 2V) + F + 2HF 



8ttGF 

For the equation of motion, eq. (128) gives 

ft + 3Hft - l -g IJ (fj - 2Vj) + Y^ft 1 + V 1 = 0. 
The null-geodesic equations are presented in eq. (156). The Boltzmann equation in cq. (171) gives 



(191) 



(192) 



\/q 2 + m 2 a- 



-cm, 



and from eq. (174) we have 



M (0 = 1 f v ^V 2 (f ffi / (c) = 1 A g 4 ^^ g /. 



(193) 



(194) 



C. Decomposed equations 

We summarize a complete set of equations necessary to analyse each perturbation type. We decompose the pertur- 
bation variables according to equation (175). Algebraic manipulations are made which can be recognized by examining 
the right hand side of the following equations. 

The scalar-type perturbation: 

A 

K-3Ha + 3ip+ — x = N , (195) 



AirGSfi + Hk+ A + ? K <p = -N u (196) 
a 2 4 



k + A+ J K X - 12kG(h + p)av = ^A- 1 V a N 2a = (197) 

k + 2Hk - AttG (Sfj, + 3Sp) + (^3H +^ja = N 3 , (198) 

X + Hx-^-a- 8irGn = ^a 2 (A + 3K ) _1 A^VV/jJV^ = (199) 

5(i + 3H ((5/i + 5p) - {fj, + p)(k- 3Ha + ^Avj = N 5 , (200) 

- i„ - > U P+ \^Lu) - -4-A- VJVi. . Nit\ (201, 
a (P+P) a a{n + p) V 3 a 2 / ^u + p 

<*A(») + 3 # (<*M(») + <*£(»)) - (M(i) + - 3i?a+ ^ Aw M^ + ~<%)o = ^5(i)j (202) 

[a 4 (M( t) + P(Q )"(«)]" 1 1 f x , 2 A + 3X TT 

IT , r a ~ "7 i \ Sp (i) + q 2 — n M ~ SI (i) 

a (/•»(») +P(i)) a a (M(i) + P(»)) V 3 a 

= j A-^Nwv = JV« (203) 

(50 + 3ff<ty - ^(50 + y 00 (50 - (k + d) - (20 + 3ity) a = JV - 0iV o , (204) 
+3i?(50 (4 ) - ^0(i) +5Z^ w ^ w ^(fe) ~ 4>(i) (K + d) - (20 (i ) + 3H<fi(i^ a = - 4> {i) N Q . (205) 

The vector-type perturbation: 
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^t^W + SttG^ +p) V M = \ (N 2a - V^A^A^) = JV&\ (206) 

*(v) + 2H y(v) _ SnGH^ = 2a (A + 2K)- 1 (v p N^ - V« A" 1 V V fj N^ ) = N { 4 V J , (207) 

feVfrft + a + » _ i (Wao _ v „ a - v , Ww) . w < : > (208) 

[o 4 (a*w+P(o )«(,>]' A + 2^ ng a i 



The tensor-type perturbation 



^7 (JV 6( i)a - V^A" 1 ^^) ee W$ )a . (209) 



-2V (Q (A + 2KY 1 (V^ 4/3)7 - V^A-^TViJV^) = iV^. (210) 
In order to derive eqs. (199,207,210) it is convenient to show 

^ (VaV, - A) (X + #X - V - a - SnGH) + 1 (a 2 *^) ' - ^G^ w 

+C« + 3ff6W - ^=^C$ - SttGhW = AW, (211) 

which follows from eq. (103). In our perturbative approach, the second-order perturbations are sourced by the 
quadratic combinations of all three-types of linear-order terms. 
For the scalar field, from eq. (121), we have: 

k 

SpW = J2 (0( fc) 4(fc) - j>t k) a-V^ k) 6cj> {k) ) + 5pW, 

k 

Q W = +p W )v W = _I £^ (fc) ^ (fc) + A- 1 V a Q(?) ) 

a k 

q(jM) = frW = Q W) _ VaA-^^Q^. (212) 

The anisotropic pressure follows from eqs. (121,177). 

In the generalized gravity theory in §IVD, eq. (128) gives 

50 1 + ZHScj) 1 - ^Scp 1 - 2a4> 1 + (j) 1 ( 30 - a - 6Ha + ^X I 
a z \ a / 

-\g IJ [FjSR + (f, L j - 2V LJ ) 6<f> L ] - ±g L \ L 5<j> L (fj - 2Vj) 

+2T JK (<j> J S4> K - A<j> J <j> K ) + Vj^J^i* + ST 1 = N g . (213) 

Equations (195-201,206-208,210) remain valid even in the generalized gravity by replacing the fluid quantities to the 
effective ones. The decomposed effective fluid quantities follow from eqs. (130,177) as: 



8ttGF 



+ 913 (¥H J - a^0 J ) + \ai.LL^ L ¥4> J - \ (f,L - 2V,l) 5<j> 1 



-3HSF + Qi? + ^jSF+ (&Ha -^ x -ZtpjF- SnG^SF 



S/j, 



(off,,) 
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8irGF 
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8ttGF 
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(214) 



where the quadratic parts follow from eq. (177). As an example, from eqs. (210,214), the gravitational wave equation 
in generalized gravity becomes 



For the electromagnetic field we can decompose: 



(215) 



(216) 



The decomposed forms of fluid quantities can be read from eqs. (149,177). Similarly, for the null-geodesic equations 
we decompose 



8e a = 5e, a + 5e^; 5e^ a = 0. 
For the temperature anisotropy, eq. (163) gives 
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r° / i \ 

*We a e% - [° C { $e a e f3 dy+ f N v dy + N T . 
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(218) 



To the linear-order this result was first presented by Sachs and Wolfe [17]; for further analyses using our notation, see 
[30]. 

For the Boltzmann equation, eq. (171) becomes 



f' + 



2 2 
~t n 



7 



,gf _ r (3)a ^7 — ^1 



* + ^^7>a + ( \a„ + »«, + C#) 7 



g^- = N c . 

aq 



(219) 



The fluid quantities can be read from eqs. (174,177). 

We emphasize that all the equations up to this point are presented without fixing the gauge conditions. In order 
to solve the equations in a given situation, we can choose any allowed gauge conditions suitable for the situation. In 
this sense, the equations are presented in a gauge-ready form. 



VI. GAUGE ISSUE 
A. Gauge transformation 

We consider the following transformation between two coordinates x a and x a 
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x a = x a + i a (x e ) = x a +c + le, b ( b . 



(220) 



The variables £ a and £ a are perturbed order quantities. To the second-order we may have 

|o = |(l)o + |(2) 0) (221) 

and similarly for £ a . For any tensor quantity we use the tensor transformation property between x a and x a spacetimes. 



~ dx^ ~ 0x c Ox^ " 

4>{x e ) = 4>{x e ), i a {x e ) = — v b (x e ), t ab (x e ) = —— 5 t cd (x e ). 



(222) 



Comparing the tensor quantities at the same spacetime point, x a , we can derive the gauge transformation property 
of the tensor quantity. We can show that a tensor quantity t transforms as [31] 



t(x c ) = t(x c ) - £ ( t + -£?t, 
where £^ is a Lie derivative along ( a . We have 



1 7 



tab(x e ) = t ab {x e ) - 2t c(a C, b) ~ U,c£ C 



1- 



+ 2t c ( a £, , b )C 4 + *cd^ C ,aC ,b + £ I 2 £ C ,(6*a)c,d + 2 *c(a£ C ,b)<j + ^afc,cd£ C + t ab ^° d I . 



Wc define 



(223) 

(224) 
(225) 

(226) 
(227) 



where £ Q is based on In terms of ( a wc set £° = C° and Q a = ( a where £ a is based on g^. Thus, we have 



(228) 



From the gauge transformation property of g ab and the definitions of our perturbation variables we can derive: 



a = a - ( e + y j - A'e - 2A u°' + A«r - B a ? 



a 2 V a a z 



(229) 



k = B a - e, a 2Ae, a - (B' a + e - B a e - B a . f3 e - B P e ia + 2c aP e 

W/3 - W/3 - —4 9 af3 ~ 2^/3,7? ~ . 



(230) 



* 7 (c^ ,/3) 



+s (Q £%) - (c' afl + 2 a -c a[ ^ e - c a f,„c - 2c 1(a c >0) + e ia e, m l?, a ?,e + 4sre°,. 



„(3)<r7 , £„(3)t7 t« + fM_„(3) ^7 , <r7 „(3) , £ (3) t 7 , (3) ,7 



(231) 
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From the gauge transformation property of T a b and using the definitions of our perturbed fluid variables in the 
normal-frame we can derive: 



1 



5p = s P - ( P > + s P ') e 6 P , a ?> + P > (ee + e, a e) + \v"ee + 1 [*Q a + & + #1 e, a , 

Q a = Qa + (Ji + P)e, a - Qfi^, a - QaA" ~ (o'a + €° + (M + P)^°, Q + + - (m' + p') £°] 

fta/3 = n Q/3 - 2n 7( «e 7 ^ ) - (n' a/3 + 2^n Q/3 ) e° - n a0 , 7 c 

+ [2Q {a + (n + p)e^ a ] e%) - [ 2 Q 7 + (m + p)e ' 7 ] c°, r 



(232) 
(233) 

) 

(234) 
(235) 



Under the gauge transformation the individual fluid quantities Sfj,^, Sp^, Q(i) a , ^(i)a(3, and S<j>^ transform just like 
the corresponding collective fluid quantities in eqs. (232-239) with all the fluid quantities changed into the ones for 
the individual one. Using the vector nature of J(j) a we have 

6i {i)0 = 51 {i)0 - (/ W o£°)' - (%)o£°)' - 6I(i)o, a e + V + l[i)of'f + [l(i)o + £\0] ' , (236) 



+1, 



(i)0 



(€°,-O , + (€ ^), a 



(237) 



The fluid quantities we use in this work are based on the normal- frame four- vector where n a = 0, see eq. (83). 
It is convenient to have the gauge-transformation properties of the fluid quantities in the energy-frame where we set 
Q a = 0. These can be derived cither by applying the frame-transformation rule presented in eqs. (87,88) or directly 
from the gauge transformation property of the energy-momentum tensor in eq. (82) with Q a = in the energy-frame. 
We have 



SfL E 

V E - B n 



5 P e - ( P ' + 5p Ei ) e s P E j a + \p n ee + P ' (ee + cej , 

v E B a +e, a (v e -B a )'e+ a - (v e B a ) e° - (v e Bp) e, a (v e - Ba ) e 

+ (v Ef} + e) (flS, 7 r + 24 3 ( U 7 ,«) + ( A ( 2 ^° - - &) + B « + 4 e °) 2C ^ e 



ir 



a/3 



n 



a/3 



n 



El 

a/3 



2 a -U 

a 



£° _ nf/3, 7 £ 7 - 2n 7 J (Q £ 7 (3) . 



(238) 



From the gauge transformation of <j> we have 



(239) 



Using the gauge-transformation property of a vector quantity k a similar to eqs. (222,225), and using the definition 
of k a in eq. (150) we can derive 



v v a 



— + — 

V V 



5v , 



_^0" _ 



a 
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I e «" _ 2 -e«' - f— - 3^ I ' 
a a z 



Using the scalar nature of the temperature T and eq. (224) we can show 

Sf(x e ) = ST(x e ) - (T' + ST') f - ST a C + T' + + \t"£°Z . 

Using the vector nature of the electric and magnetic vectors and eq. (225) we can show 

E a (x*) = E a (x e ) - E' a e - E a3 e - Ep%P , 



(241) 



(242) 



(243) 



and similarly for H a . Thus, E a and H a are gauge- invariant to the linear-order. 

Since p a = dx a /d\, under the gauge transformation we have p = p a + £ a h p b . Using the definitions of q and 7" in 
Eq. (169) we can derive 

4 = Ji + -e + ^ 2+m2a2 e al « + Vg2+m2fl V Jr ( a - e + -e) - 

la q ' q L V a / 



q 2 + m 2 a 2 



+ 



I a" :\q~ I 2/,/ V a' 2 \, n ,. 



2 a q 2 + m 2 a 2 a 2 



1 m 2 a 2 

2 g 2 



(244) 



(245) 



As 7" always appears together with perturbed order terms multiplied, it is evaluated only to the linear-order. From 
the scalar nature of / we have 



f(x e , q, 7 e ) - f(x e , q, 7 e ) = /> e + | e , q + Sq, 7 e + S^). 
At the same momentum space and the spacetime point, we have 



(246) 



sf = sf - re - fj q + ree + \ree + v[ q e^ + + - s f*c - s/^s-r. (24?) 

Using cqs. (244,245,229) we have 



si = s f - re - 9 /„ a -e + Vq2+ q m2a2 e, a r 



-of A Vg2+ g ro2 °V , a [7« (a - e) - ^ + 



A'e - A, a e 



+? ^ - a + a ? 'V + ? ? + 2 a q 2 + m 2 a 2 U a U 2 5 < Q? 
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1 V + mV o g 2q 2 + m 2 a 2 a'\ 0t0 

2 o 2 ? < Q? ^ 7 7 a 2 a 2 C C 



4°4°. Q 7 a 



yjq 2 + m 2 a 2 a 

+ 2 9 ^ |^ ^ +2 g C - 7 + ^ € ^ ^ 7 7 J 



-s.fA a - 6f, 



yjq 2 + m 2 a 2 , 



(4 -4 , P T1 ) + \€ , ~ 2? |/3 - 2^ J 7^ 



(248) 



Notice that with our phase space variables introduced in Eq. (169) the distribution function / is spatially gauge- 
invariant to the linear-order. We can check that the gauge transformation property of Sf in Eq. (248) is consistent 
with the gauge transformation properties of the fluid quantities identified in Eq. (174). 
We further decompose £ Q (and similarly for ( a ) into the scalar- and vector- types as 



1 



4a — (2^ ,a 7 ^ a ' 



(249) 



with £^ a i Q = 0. In order to fix the gauge we can impose three conditions on three variables such that these conditions 

can fix 4°, £ and £ a "\ We call these conditions fixing £°, £ and ffl as the temporal, the spatial and the rotational 
gauge fixing conditions, respectively. 

The decomposed variables in eq. (175) and others transform as 

a = a-- (a£°)' + A*, 

7 = 7- k+\ (a + ^ (3) ) (SA^V-V 3 ^ - C£) , 

= V ^4° + - ( A + 1^3))" 1 (SA-^V^ - C£) , 



£(„) = _ + 2 ("a + ii?( 3 )) (V"Q a/3 - VaA-^TV^C^) , 



c*S = - - ^ € >S - ^ (v Q V/3 - |sSa) (a + ^ (3) ) " (3A" WQ 7 , - C. 



7 



-V (a (^A + ii?< 3 )) (V^Q/3 )7 - V^A-V^V 5 ^) 

Sp = Sp- p'e + Sp^ , 

v = v-e--^— A- 1 V a Q €a , 



M + p 



(«) - „(«) 



M+P 



n = n + -a 2 A + -i?( 3 ) A^vVn^, 



n<"' = n£° + 2a f a + ii?( 3 )) (Viifa/j - VaA-^vn^) 



nS = ng - n^, - § (v Q v, - i ff g> a) (a + ^ (3) ) " A-v^n^ 
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E 



(em) 



(V^n^-VflA-^V*!!^) 



0'£° + <ty £ , 



a o z 



+A- 1 

Se^ = Se a - Se. a , 
5T = 8T- T'£° + STs, 



5f£ 
V 



a 



1/3° 



EM = £?<»> + ^ - V Q A- 1 V^£ /3 j, 



(250) 



where indicates the quadratic parts of eq. (229) and similarly for other variables. For Sf we have used /' = 
which follows from eq. (171) for C[f] = to the background order. 

Using t instead of 77 (indicated as 0) as the time variable, from the definition dt = ad-q we can show 



,0 = ^(1 



-He 



(251) 



B. Linear-order 

From eq. (250) we find that the decomposed metric and matter variables transform to the linear-order as: 



a \a I a \ a 2 



J 



Sfi = 5fi- fit*, 5p = 5p-fi£ t , v = v--£ t , ft = n. 6$ = 5(j> - </>£*, 

a 

r(") — n(«) 1 n t(v) n(v) _ (-i(v) _ e(v) _ yu(v) f (v) _ (y) fr(v) _ tt(v) 

3§ = c$> ng = n« (252) 



_Z. Temporal gauge conditions 

Temporal gauge fixing condition, fixing £*, applies only to the scalar- type perturbation. To the linear-order, we can 
impose any one of the following temporal gauge conditions to be valid at any spacetime point: 



synchronous gauge : 


a = 


- 


- £*(x), 


comoving gauge : 


v = 


- 


* e = o, 


zero— shear gauge : 


X = 


- 


- e = 0, 


uniform — expansion gauge : 


K = 


- 


- e = o, 


uniform — curvature gauge : 


if = 


- 


- e = o, 


uniform — density gauge : 


S/j, = 


- 


♦ e = 0, 


uniform — pressure gauge : 


Sp = 


- 


♦ e = 0, 


uniform— field gauge : 


5<t> = 


- 


-> £* = 0. 
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Except for the synchronous gauge condition, each of the other temporal gauge fixing conditions completely removes 
the temporal gauge mode. In the multi-component situations in addition we can choose one of the following conditions 
as the proper temporal gauge condition which also removes the temporal gauge-mode completely: 

= 0, 5p(i) = 0, v (l) = 0, 5<t> {€) = 0. (254) 

All these variables which can be used to fix the temporal gauge freedom in fact do not depend on the spatial gauge 
transformation, £, thus are naturally spatially gauge- invariant. 

The followings are some examples of combinations of variables which are temporally gauge-invariant: 

1 H H 

Sfi v = 5 (j, - fiav, ip x = tp-Hx, v x = v X, <p v = <p-aHv, ip Sl p = ip -5<j> = — -5(j) v . (255) 

O, (j) (j) 

These are completely (i.e., both spatially and temporally) gauge- invariant to the linear-order. Any variable under 
any gauge condition in eqs. (253,254) (except for the synchronous gauge), has a unique equivalent gauge-invariant 
combination. For example, we have 

<p x = <p\ x =o- ( 256 ) 

Thus, (p x is the same as <p variable in the zero-shear gauge where we set \ = 0- 

All the equations in §IV and V are presented without imposing any gauge condition. The equations are arranged 
using the above variables in eqs. (253,254) which can be used in fixing the temporal gauge condition. This allows us to 
use the various temporal gauge conditions optimally depending on the situation, thus the equations are presented in a 
sort of gauge-ready manner. Usually we do not know the most suitable gauge condition a priori. In order to take the 
advantage of gauge choice in the most optimal way it is desirable to use the gauge-ready form equations presented in 
this paper. Our set of equations is arranged so that we can easily impose various fundamental gauge conditions in eqs. 
(253,254), and their suitable combinations as well. As we have so many different ways of fixing the temporal gauge 
conditions it is convenient to denote the gauge condition, or equivalently, gauge-invariant combination, we are using. 
Our notation for gauge-invariant combinations proposed in eq. (255) is convenient for such a purpose in the spirit 
of our gauge strategy [12,29]. The notation is also practically convenient for connecting solutions in different gauge 
conditions as well as tracing the associated gauge conditions easily. Compared with the notations for gauge-invariant 
variables which were introduced by Bardeen [6,7] we have 



t m = 5 v = 5n v /n, ■$>H = tp x , vf ] =kv x , pirf>=6p, pir^ = --^11, ( = (p Sl (257) 

etc.; we ignored the harmonic functions used in [6]. The perturbed curvature variable in the comoving gauge 1Z often 
used in the literature is the same as our ip v which is the same as ips^ in the scalar field. 



2. Spatial gauge conditions 

The spatial gauge transformations £ and £^ affect the scalar- and the vector-type perturbations, respectively. Due 
to spatial homogeneity of the background we have natural spatial gauge fixing conditions to choose [7] . We have two 
natural spatial gauge fixing conditions. From eq. (252) we can see that: 

B-gauge: = 0, flW = - £(x,i)«a, ^ } (x), (258) 
C-gauge: 7 = 0, C^=0 -> £ = 0, & = 0. (259) 

For P we have considered a situation where the temporal gauge condition already completely removed £*. We call 
the spatial gauge fixing conditions in eqs. (258,259) the B-gauge and the C-gauge, respectively [5]. These gauge 
conditions are imposed so that we have 

B — gauge : B a = 0, 

C - gauge : C af3 = V g ( ^ + eg . (260) 

Apparently, the B-gauge conditions fail to fix the spatial and the rotational gauge modes completely, thus, even after 
imposing the gauge conditions we still have the remaining gauge modes. Whereas, the C-gauge conditions successfully 
remove the gauge modes. To the linear-order, the variables \ an d introduced in eq. (178) are natural and unique 
spatially gauge-invariant combinations. Notice that in the C-gauge \ is the same as a(3, and %?( v > is the same as 
Thus, /3 and Ba^ variables in the C-gauge conditions are equivalent to the corresponding (spatially and rotationally) 
gauge-invariant combinations x/ a an d \ respectively. 
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C. Second-order 



1. Gauge conditions 

If we use any of the gauge conditions which completely fixes both the temporal and the spatial gauge modes to the 
linear-order, the gauge transformation properties of the second-order variables, say ip^ in cq. (183), follow exactly 
the same forms as their linear counterparts. Using the transformation of 8<j) in eq. (239) as an example, to the linear 
order we have 

5$W =6<f>W -<^ 0(1) - (261) 

If we take gauge conditions which remove (fix) £° and £ Q completely to the linear order we have C 0(1) = = 
Thus, from cq. (239) we have 

5$^=S<f>^-4>'^ 2 \ (262) 

which shows exactly the same form as in eq. (261). Thus, the gauge conditions in eqs. (252,254) apply to second- 
order perturbation variables as well, and we can impose similar gauge conditions even to the second-order. For 
example, in the zero-shear gauge we impose x = as the gauge condition to the second-order, thus x^ = = X^> 
otherwise mentioned, we always take the C-gauge for the spatial and rotational ones. In this gauge condition the 
gauge transforamtion properties are completely fixed, and the gauge-modes do not appear. Thus, we anticipate that 
each variable in that gauge condition has unique corresponding gauge-invariant combination of variables. Thus, using 
ip, we have </5| x =o, c-gauge is free of the gauge-modes. We denote the corresponding gauge-invariant combination as 

Vx (263) 

with the C-gauge condition assumed always. To the linear-order we have <p x = <p — H\, but to the second-order we 
need correction terms to make ip x gauge-invariant. Construction of such a gauge-invariant combination will be shown 
below. 

From eqs. (57,58), assuming pure scalar mode in the C-gauge we can show 

K a p = -{l-a) x M0 + 2x, {a <p,p) ~ \g ( $ [- (1 - «) Ax + 2x'V 7 ] , (264) 



/"J z 



(265) 



Thus, the gauge condition x = implies K a p — 0, justifying its name as the zero-shear gauge to the second-order. 
Similarly, the gauge condition ip = implies = -^R^ (we also have R^ — g-R^S^g), justifying its name as the 
uniform-curvature gauge to the second-order. We can show that the names of gauge conditions in eq. (253) remain 
valid to the second-order. 

In the perturbative approach, apparently, this method can be similarly applied to any higher-order perturbations. 
As long as we work in any of these gauge conditions, the gauge modes are completely removed and the behavior of 
all the variables is equivalently gauge- invariant. As the variables are free of gauge mode, these can be considered as 
physically important ones in that particular gauge conditions we choose. We can also choose different gauge conditions 
in the second-order compared with the ones imposed to the linear-order. Examples will be shown below. 



2. Constructing gauge-invariant combinations 

Let us explain a method to derive the gauge-invariant combinations using an example. Since the gauge transfor- 
mation properties of 5fj, and 5<p are available in a convenient forms in eqs. (232,239) we consider the gauge-invariant 
combinations involving these two variables to the second-order. Thus, we consider the case with a scalar field. To the 
linear-order we can construct various gauge invariant combinations involving 5fi, and as examples we consider two 
cases 

Sfj,s<p = 5 (j, - T7<$0, <W = - (266) 
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Apparently, the combinations in eqs. (266) are not gauge-invariant to the second-order. In order to construct the 
gauge-invariant combination in the gauge with Scf)^ = we construct 5fi — ^S<t> using eqs. (232,239). Then, in the 

RHS we have quadratic combination of linear-order terms involving £° and £ a . As the spatial and rotational gauge 
we consider the C-gauge conditions which remove the corresponding gauge modes completely. This can be achieved 
by taking 

= - (V + C* (t,)a ) +r a + C (v)a , (267) 

which follows from eqs. (252,249), and moving terms with overhat to the LHS. Now, coming to the temporal gauge 
freedom, if we want to consider the uniform-field gauge we take 

e = ~(n-H), (268) 

which follows from eq. (252), and move terms with overhat to the LHS. Then we have a gauge-invariant combination 



,2' 



= ^-^ + ^g. (269) 

We have (J/i^ = "^A* | ,5^(1) — 0— 50( 2 ) c-gaugc' ^ nus ^A 4 *^ i s the same as <5/x under the gauge conditions 5(j>^ = = 5<j)^ 
and the C-gauges. If we want to take the uniform-curvature gauge to the linear-order we take 

? = ~{<P-<P) (270) 

which follows from eq. (252), and move terms with overhat to the LHS. Then we can identify a gauge-invariant 
combination 



This combination is equivalent to 5/j, in the following gauge conditions: 5<fi = in the linear and the pure second-order 
part [i.e., d/i^ — (/x' / ' 4>')8(j)^], and = in the quadratic parts, and the C-gauges. By replacing the linear-order 
part of eq. (271) with - </>'(a/a')<yS (1) we can make another gauge-invariant combination 

*/•„...,»» = - *'^if m + «(' (31 - - («f - (r° + C W °) - J (if- £<V)'v 

- (£)' 7 (* - 5*?») * - 4* 7 (**" 3^)" (272 » 

which is the same as Sfj,\ v (i)_ 0= ^(2) c-gauge' ^ ne calculation becomes simpler if we take the C-gauge condition: this 

sets 7 = = Ca"^ , thus we can simply set £ Q = (£ = = ) . Similarly we can construct diverse combinations of 
the gauge-invariant variables: several useful gauge-invariant combinations will be presented in the next subsection. 

In the following, as in eq. (269), a gauge-invariant notation, say, ip v indicates a combination which is equivalent 
to if in the comoving gauge (v = 0) to all orders, thus, = = v^ 2 \ and in the C-gauge. In order to denote 
gauge-invariant combinations valid to the second-order, we introduce the following notation 

<p v = <p - aHv + ipM , Vx = ^-H x + ^\ 5fjL v = 5ii- iiav + 8^\ v x = v-- X + v x q \ (273) 

a 

etc., where the upper (q) index indicates the quadratic combinations of the linear-order terms. In the following we 
always take the spatial C-gauge. We note that 6[i^ is the quadratic correction terms to make 5^ v a gauge-invariant 
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combination to the second-order; thus it differs from, say, S/i^ — fiav^ q \ As ip v is the same as <p in the v = gauge, 
we have ipi q ^ vanishes under the v = gauge, i.e., 

fv I" ~~ fv \v=0 — u - 

Using the definition of our gauge-invariant combinations we can show, for example, 



ip v = ip- aHv + ipW = ip x - aHv x + ip^ 



X' 



where in the second step we have evaluated the first step in the zero-shear gauge. Thus 

<pW \ x = Vv -{!P x - aHv x ) = <pM - aHv^) , 

and similarly for other correction terms. 



(274) 
(275) 

(276) 



3. Gauge-invariant variables 

Now, we present several useful gauge-invariant combinations explicitly. We assume = and pure scalar-type 
perturbations. We take the spatial C-gauge, 7 = 0. As long as we take the temporal gauge which fixes £° completely, 
we can set £ a = 0. The metric becomes 



A = a, B a = -x, a , C a f3 = (pg ( 3- 



From eqs. (229-235,239) we have 



a = 



a \ a J 2 



0/ 



2 \ a a 



a 



a or 



\ - ¥°- a ^) - W 1 ™' - ¥'-) 



--A -1 
2 



2 \a 

X" 

k = k + I 3— — 6— H £ + quadratic terms, 



-2< a - - f x ' + -x) e° - -x, Q e 0/ + ee 



a ? < Q 2 ? ? 



a J a 



1 



-A V c 



2 \ a /i + p 

a A* + P a M + P 

1 



, « A* + P 
4 h 



a /j, + p 



(277) 



(278) 



In the transformation of v we have used eq. (201). We have ignored writing the quadratic terms in the transformation 
of k; this can be read from the definition of k 



1 



K = 

a 



2a 



a 



(279) 



which follows from eqs. (179,57). 
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Following the prescription in the previous subsection, from eq. (278) we can construct the following gauge-invariant 
combinations 



fx = <P- 



H X {<p x + 2H lfix ) X - \ (H + H 2 ) x 2 + ^ [r a X., a - A-'V a V^ (xm)] 



+HA- 1 \7 a [2a xX , a + ( X ~ H X ) xA + ^HA' 1 



^X^a-^A^WV ( X ,aX,t 



tp v = tp- aHv - (tp v + 2Htp v ) av- - I H + 2H Z - H 



P-+P\ 2 2, 1 



1 

la" 



A _1 V Q V /3 



(2x - av) a v t0 



+ aHA- l S7 a 



H + p 
H + p 



a v + — (2y — av)' v , 
4a 



-v. 



1 1 

a 2 [i + p 



A Sfl v lflH 2 2 /U+Pa 

o„ = o aw an -aw w u, 

^ jU 2 Ji jU 



+ ^aA- 1 V« 



1 1 



-v., 



1 

w x = w - -X- 



fi + p a* n + p 
fi+p 



(280) 



(281) 



(282) 



+ \H + 



fi + p 



X + 



2a 



+ ^A- 1 V a 
a 

+ 4 A 



fi + p 

2a x x, a + X,a (X - H X ) - -P-TZX 



fi + p" a 2 fi+p 
1 x - x ^ a - \A-^ a ^ { XM ) 



(n^-ofAn) x ^ 



(283) 



Xv = X- aw - ^a 2 \ 11 
+aA- 1 V a 



M + P \ „,2 
fi + p 



-2a v v tCt - {Xv + H Xv ) a v- X, a v - 



6fi v 

1 

fi + p 



1 1 

a 2 fi + p 



(n^ a -«yfAn)t, i( 



(284) 



To the linear-order, S v (equivalently, S in the comoving gauge) behaves like the Newtonian density perturbation, and 
v x and — (p x (equivalently, v and — (p in the zero-shear gauge) behave like the Newtonian velocity and the gravitational 
potential. Also to the linear-order ip v is known to be the best conserved quantity in the super-sound-horizon scale. 
For extensions of these results to the second-order, see §VIIC and VII D, respectively. 
In the case of a scalar field we have 



H 



1 



1 / a 2 H 



+ - 



2a 2 c/> 



2<f) 



2<p 



5<f> v = 54>--<p--6<l> v <p + ^ Wi 



a 2 H 



2a 2 H 2 



(285) 



(286) 



Thus, to the second-order we have ipg^ ^ — ^5(j> v . By evaluating the RHSs of eqs. (285,286) in the ip = gauge and 
5<j) = gauge, respectively, we have the following relations between the two gauge-invariant variables 

1 /a 2 /T 



<Ps<i> = — j-o<p v + -j-d(p v d<p v 



2a 2 <p 



+ 



2a 2 (f> 



2<P 



2(p 



(287) 



4> 



1 



f = —jjVW + ^WW* + [ — ) 



H 2 



2a 2 H \ H 
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2a 2 H 2 



Xs<p Vd<l>,a 



2H 



(288) 



To linear-order 8(j) v (equivalently, 8<p in the uniform-curvature gauge) most closely resembles the scalar-field equation 

in the fixed cosmological background metric [14]. Since 8<j) = implies Qofi = (thus = 0), we have ips,p — <p v ; 
from eq. (114) we notice that this is valid to the second-order. 



A covariant density gradient variable 



4- Spatial gradient variable 



A a = -h h a fiM 7 
/i 



(289) 



is gauge-invariant to the linear-order [32]. In [32] the energy-frame is taken. Using the u a -frame in eq. (53) we have 



6' + — (5 — A) 



(Kf (290) 



where we set A a = A a with A a based on grj; from u a A a = we have A = —V a A a . From eqs. (85,86) we notice 
that 5 is frame invariant to the linear-order, and we ignore superscript E in such cases. Using the prescription in eq. 
(88) we can express A a in the normal-frame as 



1 



6 N „ + 



N 



1 + 8 



H pi + p 



8<+>-{8-A) 



Qa (Q N(3 Qp).a p! 



1 



h+p h(h + p) n ip + p) 2 



[{Sn + 8p) + H a/3 Q 



(291) 



Thus, under the condition Q a = we have A Q = <5, Q /(l + <5); Q a = can be achieved by the comoving gauge condition 

(Q = 0) and the irrotational condition (Qa^ = 0); see [33] and the note added in proof of [21]. Under the gauge 
transformation, either using eqs. (232-235) for eq. (291) or using eq. (238) for eq. (290) we can show 



A« = A Q AU - (^f). a - 2A [a ,0f 







(292) 



Thus, A Q is not gauge-invariant to the second-order. To the linear-order, the scalar-type part becomes A Q = 8 Vt , 
where the gauge-invariant combination 8 V is the same as 8 in the comoving gauge. 



VII. APPLICATIONS 



A. Closed form equations 



From eqs. (196,197), eqs. (197,200,201), eqs. (199), eqs. (199,201) and eqs. (195,197,199) we can derive, respec- 
tively: 



8 fi v + 3H8fi v - A+ J K [a(p+p)v x + 2HU] = Sffl + 3H8fi v g) - A + J K a(p + p)v<?) 

+N 5 + (p+p) (A 2 (s) + 3aHN^ , 
Vx + a x + 8nGU = <pM + - jvj s) , or <p x + a x + 8ttGI1 x = - jvg> , 

x x a\ M + P 3 a 2 p + pj x x a \ x p + p 

if x + H<p x + 4irG(p+p)av x + 8wGHIl = + Hp® + AixG{p + p)av<?> + | (n - - HN^. (297) 



N, 



6 ' 



(293) 

(294) 
(295) 

(296) 
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These equations are presented using mixed gauge-invariant variables. We note that to the linear-order, Sfx v , — ip x , 
and kv x {^ Vv x ) closely resemble the Newtonian density perturbation, perturbed gravitational potential, and the 
perturbed velocity perturbation, respectively [34,35,13] . To the linear-order these equations were presented by Bardeen 
in 1980 [6]; see eqs. (4.3,4.8,4.4,4.5,4.7) in [6], and compared with our notation, see eq. (257). Using eq. (297) and 
eqs. (293,296,297) we can show 



$ = fv - 



K/a 2 



4ttGV + P) 
H 2 



2H Z 



n 



AttG(/i + p)a 
4ttG(m + p)a 2 V x ^ x ) ft + p \ 



fi + p 
2 A 



3 a 



+ + 
n | + <5> {q) + N, 



<i>' 



(298) 
(299) 



where 
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We have introduced an entropic perturbation e by 



5p=c s 8n + e, c s 



(300) 



(301) 



Defined in this way, e is not necessarily gauge-invariant to the second-order. To the linear order we have e = tt 
introduced above eq. (24). Combining eqs. (298,299) we can derive 
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where we used 
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(304) 



[i> in eqs. (302-304) differs from the perturbed velocity related variable used in the rest of this paper.] The equation 
using v in the linear theory was first derived by Field and Shepley in 1968 [36], see also [37,21]. Using eq. (293), eq. 
(303) gives an equation for S v . Using eqs. (293-296) we can derive equation for 5 V in another form 
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The above set of equations is valid for a general imperfect fluid. A minimally coupled scalar field can be regarded 
as an imperfect fluid with special fluid quantities. We additionally have an equation of motion of the field which is 
actually included in the energy and momentum conservation equations. In fact, the above set of equations is valid 
even for multi-component fluids and fields. In such cases, the fluid quantities become collective fluid quantities and we 
additionally need the energy and momentum conservation equations for individual fluid and the equations of motion 
for the individual field. 

In the single scalar field case, from eqs. (212,190,293,301) we have 



1 - c 2 s A + 3K 
4nG a 2 



'tp x -tpM)+N e 



N e ee L-^ QjV! - HN^ + Sp^ - + 3H(l - C 2)aA" 1 V a gi 9) . (306) 



Equation (298) remains valid, and eq. (299) becomes 

Hc 2 A A 



$ = 



- 2 - -^- p ( Ne + + ^ + N * (30?) 



4t:G((i +p)c 
where 

c\A ee A + 3(1 - c 2 s )K. (308) 

Therefore, eqs. (302,303,304) remain valid with c s and e replaced by ca and N e ; in eq. (302) one can show that we 
can ignore the operator nature of A -1 in the c\. 

The rotational perturbation and the gravitational wave are described by eqs. (208,210), respectively 

rfr a+2A ' n£ " + <, (309) 



a 4 (fi+p) 2a 2 fi+p 

SttGU^ + N^, (310) 



Mt) , oTfMt) A-2K (t) _ 1 



where v^p = Cq/3- We note that all the equations above in this section are valid for general K. 



B. Solutions to the linear-order 



1. Scalar-type 



We consider a single component ideal fluid. Several known solutions in the literature are the followings. 
(i) In the large-scale limit (the super-sound-horizon scale), i.e., ignoring c 2 A term compared with z"/z and 
(l/z)"/(l/z) terms in eqs. (302,303), we have general solutions 1 

mt) = cik)- d{ k)^f^^ )d t, (3ii) 

<p x {k,t) = 4irGC(k)^- J a ^ 2 P ^ dt + d(k)^. (312) 

C(k) and d(k) are the coefficients of the growing and decaying solutions (in an expanding medium), respectively. To 
the second-order in the large-scale expansion we have 



1 k is the wave vector with k = |k|. With the wave number k appearing in the equation the variables can be regarded as the 
Fourier transformed ones. To the linear-order each Fourier mode decouple from the other modes and evolves independently. 
The same equations in configuration space remain valid in Fourier space as well. Thus, we ignore specific symbols distinguishing 
the variables in the two spaces. Only in this subsection concerning the linear theory we use the Fourier transformation. 
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k 2 n drj 



m r 



^ x = 4ttGC— I z z dn + d— \ l + k 



/r] 1 / rn \ /-J7 (-r; 



(313) 
(314) 



We emphasize that these solutions are valid for general K and A, and time-varying equation of state. 

(ii) In the small-scale limits [c 2 fc 2 3> z"/z, (l/z)"/(l/z)], if we further assume that c s is constant in time, eqs. 

(302,303) give the general solutions: 



v = z$ cx e ±iCskv , u 



1 



ip x cx e 



±ic s kr) 



(315) 



(iii) For K = = A and w = p/fi = constant we have exact solutions [6,38]. For the background, from eqs. (184,186), 
we have 



a cx i 3 <i+™) cx 771+3™ , aHrj 



1 + 3w 

Equation (304) gives z cx z cx a, thus 

z" _ 2(1 -3to) 1 (1/z)" _ 6(1 + w) 1 



z (l + 3w) 2 r7 2 ' (1/z) (l + 3w) 2 r7 2 ' 
In this case, eqs. (302,303) become Bessel's equations with solutions: 

3(1 -w) 



v = z$ oc s /r}(J v (x),Y u (x)) , x = c s k\ri\, v 
1 



:ip x oc ^/rj(Jp(x),Y p (x)) , v 



2(1 + 3w) 
5 + 3w 



We have v = v + 1. Using eqs. (298,299) we can normalize the solutions as 
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Equation (293) gives 
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In the large-scale limit (a; <C 1) we have 
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(319) 
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(322) 
(323) 



where for c = 0we have additional 2 In a; factor in the C2-mode. By matching with the general large-scale solution in 
cq. (311) we can identify 



ci = 2T(z/+l)C, c 2 = - 



1 



3(1 + tu) 2T(i/) a 2 r7 



d. 



In the large-scale limit (1 < 1) we have 

$ cx C, da-^ 1 -^, 

5+3v> 

ip x cx C, da 2 , 
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Equation (327) follows from eq. (293) which gives S v cx a 1+3w ip x cx rj 2 ip x in general. Equation (327) includes the well 
known solutions in the matter (w = 0) and radiation (w = eras [1]: 

mde : S v cx Ca, da~i cx Ct% , di -1 cx Cry 2 , d?7~ 3 , 

rde: 5 V cx Ca 2 ,da^ x cx Ct, dt~i cx Cn 2 , drf 1 . (328) 
If we consider only the C-mode which is the relatively growing-mode in an expanding phase, we have 

$(x,t)=C(x), (329) 
^(x,t) = !±|£c(x). (330) 

The nontransient mode of $ remains constant in the super-sound-horizon scale. Whereas, the one of ip x jumps as the 
background equation of state changes. Still, it is — <p x , not — ip vi which closely resembles the perturbed Newtonian 
gravitational potential [13]. 

The asymptotic solutions (i) and (ii) remain valid for the scalar field with K = 0; in this case we have c 2 replaced 
by 1. The background solutions for w = constant case considered in (iii), eq. (316), are valid for a scalar field with 
exponential potential of the following form [39] 



V ~ l2*G(l + wy e ' 6^G(1 + W ) lnt - [661) 

For the perturbation, from eq. (308) and the prescription below it, the same equations of the fluid remain valid with 
the coefficient of Laplacian term replaced by 1 (for the field) instead of c 2 (for the fluid) [14]. Thus, the perturbation 
solutions in the fluid remain valid for such the scalar field with x = k\r}\, instead of x = c s k\i]\ for the fluid case. 

We emphasize that, if a solution is known in a given gauge condition the rest of the variables in all gauge conditions 
can be derived through simple algebra; such solutions are presented in tabular forms for an ideal fluid and a scalar 
field in [38,40,15]. Solutions in the situations of generalized gravity theories considered in §IVD can be found in [24]. 
The cases with multiple components of the fluids and fields were analysed in [41,42]. 



2. Vector-type 

The rotational perturbation is described by eq. (309). If we assume 

n£° - o, we have a general solution 

a. a 3 ( / i + p).^)(x,t)=i(:)(x). (332) 

Thus, to the linear-order the rotational perturbation is described by the conservation of angular momentum, and is 
transient in expanding media. We note that eq. (309) follows from the T^. b = 0, thus independent of the gravitational 
field equation. Thus the presence of scalar fields or the generalized gravity theories considered in §IV D do not affect 
the vector- type perturbation of the fluids [23] . 



3. Tensor-type 

Now, we consider the gravitational wave with K = and IT^ = 0. The basic equation is presented in eq. (310). 

(i) The general large-scale (k 2 <C a" /a) solution is 

C ( :l{k,t) = c%{k)+d%{k) f^. (333) 

Thus, ignoring the transient mode (d a p) in an expanding phase the tensor-type perturbation is characterized by its 
conserved amplitude c^'l(fc). 

(ii) In the small-scale limit (k 2 3> a" /a) we have a general solution 

c2(M)cx (334) 



51 



Thus the gravitational wave redshifts away. 

(iii) For K = = A and w = constant, we have an exact solution: 



Thus, we set 



v$, = aC§cxJrj{J v {x),Y v {x)), x=k\ V \. (335) 



036) 



In the large-scale limit (x <C 1) we have 
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iitiona 

in eq. (333) we can identify 
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where for v = we have additional 2 In a; factor in the c^-rriode. By matching with the general large-scale solution 

c (t) -2T(v + l)c {t) c (t) 7T rjx 2v (t) 



Thus, 

- cS,4>-K^) cx cS,4>"^ K cS,4V^- (339) 



C. Pressureless irrotational fluid 



The equation of S v was derived in eq. (305) or can be derived from eqs. (303,296). In the pressureless case, a 
simpler route is to use the basic equations in eqs. (195-201). 

We consider a pressureless fluid, thus Sp = = IT Q( g , and ignore the vector-type perturbation. For the spatial gauge 
we take 7 = 0, thus j3 — x/ a - If we take the temporal comoving gauge (v = 0) we have Q a = 0. Equation (201) 

(s) 1 

gives a v = —aN^J = — ^f3 v ,aP v ' a , thus a vanishes to the linear-order; in the pressureless medium, to the linear-order, 
the comoving particle follows geodesic, thus v = implies a = 0. From eqs. (200,198), first evaluating these in the 
comoving gauge, we can derive 

S v = k v + -V • (S v Vv x ) , (340) 

k v + 2Hk v 4ttGhS v = 1 [(Vv x ) ■ (Vv x y a } ]a + C§ (d^ + JrX/^) , (341) 

where we have used \ v = —av x following from eq. (255), and k v = ^Aw x following from eq. (197) with K — 0, both 
to the linear-order. 

In order to compare with Newtonian analysis we introduce u = — \7v x to the linear-order. By combining eqs. 
(340,341) we have 

S v + 2H6 V - 4irGfl5 v = ~ [oV • (<5„u)]' + • (u • Vu) + G% |V«^ + C {t)afi ^j . (342) 

We note that, to the linear-order, the growing solution of the gravitational wave remains constant in time in the super- 
horizon scale, whereas it redshifts away (C^l oc a -1 ) in the sub-horizon scale; see §VIIB3. Ignoring the gravitational 
wave we reproduce correctly the corresponding Newtonian equation 

S + 2HS- 4nGgS = ~ [oV ■ (Su)\ + V • (u • Vu) . (343) 

We note that our eq. (342) is valid in the super-sound-horizon (Jeans scale) which is negligible in the pressureless 
medium, thus valid even in the super-(visual)-horizon. In the Newtonian context eq. (343) is valid to all orders in 
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perturbation, and follows from the mass conservation, the momentum conservation, and the Poissons' equation given 
as [10] 



6 + -V- u = — V • (Su) , 



a a 



(344) 



u + Hu + -V5$ = — u • Vu 



a a 



(345) 



4v 2 <5$ = 4nG§5. 



a 



(346) 



To the linear-order these equations can be compared with the relativistic version in eqs. (294,296,293) with eq. (295). 
To the second-order, however, although the final result in eq. (342) coincides with the Newtonian one in cq. (343), 
we notice some difference between eqs. (340,341,294-296) with eq. (295) and eqs. (344-346). From eq. (197), to the 
second-order, we have 



Since the RHSs of this equation and of eq. (293) do not vanish we cannot directly relate — <p x and — Vw x (or <jA~ 1 Vk v ) 
to the Newtonian counterparts 5$ and u, respectively. Still, we emphasize that the final equation identified in eq. 
(342) coincides exactly with the Newtonian one in eq. (343). For a similar conlusion in the relativistic situation, see 
[43]. For analyses of eq. (343) in the Newtonian context, see [44]. 

Using (i = as the spatial gauge condition Kasai [45] has derived a different equation compared with ours in eq. 
(342). Kasai [45] took both the comoving gauge v = and the original synchronous gauge which takes a — — j3. 
As we have shown above eq. (340) in a pressureless medium, the comoving gauge v = implies a = — i/3 Q /3' Q , thus 
vanishes if we take (3 = as the spatial gauge condition. However, in that gauge condition (the spatial _B-gauge) the 
spatial gauge-mode is incompletely fixed. Thus, comparison with the Newtonian analyses is not tranparent in that 
gauge condition. 



The general equation of the pressureless and irrotational ideal fluid can be derived from eqs. (10,12,13). The 
pressureless ideal fluid implies S a p — 0. We take the temporal comoving gauge condition, v = 0. Together with the 
irrotational condition we have Q a — 0, thus J a = 0. Equation (13) gives N >a = 0; if we use the normalization in eq. 
(55), we have N — a. Equations (10,12) give 



Apparently, the spatial B-gauge condition, B a = 0, leads to N a = 0, thus simplifying the equations. However, such a 
gauge condition leaves the spatial gauge-mode removed incompletely. We prefer to take the spatial C-gauge condition 
which fixes the spatial gauge modes completely, in this way the analyses can be equivalently considered as spatially 
gauge-invariant ones. From eqs. (348,349) we can derive eq. (342) to the second-order. We notice that, in contrast 
with the Newtonian case, in general we anticipate to have infinite perturbation series expansion, and eq. (342) looks 
valid only to the second-order. If we have the higher-order terms nonvanishing, these can be regarded as purely 
relativistic effects. 




(347) 



1. Nonlinear equation based on 3 + 1 formulation 




(348) 



(349) 



2. Nonlinear equation based on 1 + 3 formulation 



Assuming pressureless condition, eqs. (26,27,28) in the energy-frame (q a = 0) become: 



a a = 0, 

6 + \e 2 + iirGfi - A + 2 (, 



(350) 
(351) 



~2 - 2 

a - lu 



)=0. 



(352) 
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In the energy- frame the frame- vector follows the possible energy flux, thus, the energy flux term q a vanishes. Equations 
(350) and (352) can be combined to give 

^ -AirGfi + A-2(a 2 -u 2 ) = 0. (353) 

If we set /i = yu(l + 5), where \i is the background energy density, equation (353) becomes 



-.2 



3 n 3 f + o 



A\ _ i /A 

fij 3 



2 



- AttGh + A 



= 0. (354) 



This is a completely nonlinear equation. 

Now, we assume an irrotational fluid. In the energy-frame, the comoving gauge condition leads to u a = 0; this is 
equivalent to taking the normal-frame with vanishing energy flux. The momentum conservation equation in eq. (351) 
implies that our frame-vector follows geodesic path. In the comoving gauge eqs. (351,69) lead to A = —^B a B a to 
the second-order. Thus, to the linear-order we have A = which coincides with taking the synchronous gauge. We 
have 

l = jl, a u a = d t jl+-il, a B a - Ji = d t fJ,, 5 = d t 5+-5, a B a . (355) 
a a 

Only in the comoving gauge condition the covariant derivative along u a simplifies to the second-order as in eq. (355). 
In this derivation the pressureless condition is used essentially. Thus, in this comoving gauge equation (354) becomes 

5 + 2HS- AnG^S = \ (a 2 5S) ' + U 2 - \ {a5 ta B a )' - -5, a B a + 2ct 2 , (356) 

& V / o Qi CI 

to the second-order. From this we can derive eq. (342) where a 2 follows from eqs. (16,57,71). 

3. Gravitational wave as a source 

From eq. (342) considering the pure gravitational wave as the source for density perturbation we have 

S v + 2HS V - 4ttGhS v = C ( *lc {t)afi = S. (357) 

In the matter dominated era an exact solution is given in eq. (336) with v = |; in the large-scale limit, considering 
the relatively growing mode, S vanishes, whereas, in the small scale limit it decays proportional to l/(at) 2 oc a -5 , see 
eq. (334). If 5 g and Sd denote two linear-order solutions, the general solution can be written as 

9 J '5 g ( x ,t')S d ( x ,t')-S d ( X ,t')8 g (K,t') 

The particular solution is proportional to a~ 5 t 2 oc a~ 2 , and decays more rapidly even compared with the decaying 
mode in the linear theory which behaves as t" 1 . 

D. Pure scalar-type perturbation 

Equation (302) can be written as 

{ <<ii*! \i - - Ni + („ + "n)l } ' _ e£ U - - 2H* JL - N t ) . (359) 

In the large-scale limit, if we could ignore the second-order spatial derivative terms, we have 

&_$(g)_jV A + _g_ eoc (360) 
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Now, we consider K — and ideal fluids, thus e = 0. From eq. (298) we have $ = ip v . Ignoring the second-order 
spatial derivatives, from eqs. (300,195-201,99-105) we can show that N§ — In the comoving gauge we have 

$(«)| v = = 0, see eq. (274). Thus, we have 



<p v -tf = C(x) + d(x) f {f ?*f )a3 dt- : ™ 1 ! 



Therefore, ignoring the transient mode in the expanding phase, we have 

ip v -<p 2 v = C(x), (362) 

which remains constant even to the second-order in perturbations [46,47]. Thus, <p v is conserved to the second-order 
in the large-scale (super-sound-horizon) limit. 

Equation (359) is valid for p ^ 0. For p = 0we have a simpler form in eq. (299) which gives 



$_$(9)_jV i + -(e+ n)= 0. (363) 



H ( 2 A 

Equation (360) includes this as a case in the large-scale (super-sound-horizon) limit. Thus, the above results in eqs. 
(360-362) remain valid for general p. 

From eq. (114) we notice that for a minimally coupled scalar field 8<j> = implies v = to the second-order, thus 

fv = ¥>5<j>, (364) 

and the uniform-field gauge coincides with the comoving gauge. Thus, the above analyses are valid even for a minimally 
coupled scalar field. 



E. Pure rotation 



In the case of pure rotation eqs. (208,105) provide a complete set for a single component fluid; for the multi- 
component case see eqs. (209,107). Assuming = we have 
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(365) 



As a simple exercise, using eq. (250), one can check the gauge transformation properties of both sides. In the C-gauge 
condition {C a v ^ = 0) we have — which is gauge-invariant. From eqs. (332,206), to the linear-order we have 



\&a^ oc a 2 (fi + p)v [ a v ' a a 2 , thus from eq. (365) we have 



a 4 {n+p)v, 



(v) 



= a 



a d -a 4 (fi + p)K 



(v) 



cx a 



Thus, the additional second-order perturbation sourced by the RHS of eq. (365) behaves as 



a\n + p)vW = (x) + a 3 • a 4 ( M + p) iV^ 



(366) 



(367) 



The time dependent nonlinear solution is proportional to J dt/a 3 ; for w = constant it is proportional to a 3 ( 1 ""'Z 2 , 
thus it always decays (in expanding phase) for w < 1 . The lower bound of integration which could give a temporally 

constant nonlinear solution can be absorbed to . 

As we explained in §VI B 2, to the linear-order, the C-gauge condition removes the rotational gauge-mode completely, 
whereas the B-gauge condition fails to fix it completely That is, even after imposing the gauge condition we have 
some modes which are coordinate effect; under the -B-gauge, from eq. (252) we have = £q^(x). Then, in eq. 
(365) we notice an ironic situation where the -B-gauge condition gives vanishing quadratic terms, whereas these terms 
do not vanish in the C-gaugc condition. That is, although we anticipate the nonlinear solution in the C-gauge in eq. 
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(367) is physical, in the -B-gauge condition the RHS of eq. (366) vanishes, and we do not have the nonlinear solution 

in eq. (367). We can check this situation by using the gauge transformation property of Va^ variable in the two gauge 
conditions. 

Considering pure vector-type perturbation, from cqs. (175,177,234,249) we have 



(368) 



Now, let the variables with hat and no-hat correspond to the ones in the B- and the C-gauge conditions, respectively. 



As the £i"^s appear in quadratic combination, we need it only to the linear-order. From eq. (252), we have B^ 1 — 



(V) _ 



Ba + a£a ' . Since the hat indicates the B gauge, we have a£a ' = —B„ , thus = — j (Ba /a)db. Thus, eq. 
(368) gives 



B— gauge 



a\n+p)v^- \a 3 -a 4 (n + p)N& 



dt 



(369) 



C- gauge 



Therefore, in the i?-gauge the nonlinear solution in the C-gauge in eq. (367) disappears exactly. We note, however, 
that the solution in eq. (367) is physical (gauge-invariant) one in the C-gauge. 



F. Pure gravitational wave 

In the case of pure gravitational wave eqs. (210,103) provide a complete set. In the large-scale limit, thus ignoring 
second-order spatial derivative terms, and assuming K = and ng = 0, we have 

C (t) ; + 3iK7 (t) ; = Aft, 

Aft - A 4/3 Q - A- 1 (V Q V 7 Ay + V/jVJV) + \ (A- 1 V a V /J + 5$) A^V^A^, 

N$ = 2 [c^C^ - \s^C^i s C^j . (370) 

Notice that in this large-scale limit we have C^ = constant as the relatively growing solution (in expanding phase) 
even to the second-order perturbation. In this sense, ignoring the transient mode in expanding phase, the amplitude 
of C^l remains constant even to the second-order in perturbations [46,47]. 



G. Action formulation 



We consider the action expanded to the second-order in perturbations which will give equations of motion to the 
linear-order in perturbation [48,37,9]. We consider the action for a scalar field in eq. (108). The perturbed action can 
be derived by using eqs. (94,98) and the ADM quantities presented in §111. To the background order, ignoring the 
surface terms, we have 



Sbg = 



1 



16ttG 



-6 



6A 



+ 16ttG 



-V 



dtd 3 x. 



To the second-order perturbation, ignoring the surface terms, the pure gravitational wave part becomes 
1 



w 



16ttG 



(371) 



(372) 



This action is valid for arbitrary number of scalar fields and fluids with vanishing tensor-type anisotropic stress. Now, 
we consider the pure scalar-type perturbation. We assume K = 0. To the second-order perturbation, ignoring the 
surface terms, we have 
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(373) 



In this case we used the linear-order equations of motions, thus it is an on-shcll action. In the second step we used 
cq. (255). 

Maldacena has considered the perturbed action to the third-order in perturbations which is needed to have equations 
of motion valid to the second-order [47] . For the temporal gauge he used two gauge conditions, the uniform-field gauge 
{8(f) = 0) and the uniform-curvature gauge (<p = 0), and the C-gauge for the spatial and rotational ones. Compared 
with our notations we have 

CMaldacena = ¥><50 ~ Ca<£i '/'Maldacena = 8(j)^ . (374) 

Thus, CMaldacena is conserved in the large-scale limit, see eqs. (362,364). To the linear-order, from eq. (255) we have 

V5<t> = ~(H/(j))S(j) v , thus CMaldacena = - {H / <f>) CMaldacena- 

VIII. DISCUSSIONS 

We have presented the basic equations to investigate the second-order perturbation of the Friedmann world model. 
In order to serve as a convenient reference for future studies and applications we have presented some useful relations 
and quantities needed for the second-order perturbations. The present study is, apparently, not entirely new in this 
rich field of cosmology and the large-scale structure formation. In the late sixties Tomita has presented a series of work 
on the subject in the context of a fluid [49]. Studies in the context of the ideal fluid or the minimally coupled field 
can be found in [31,50,47]. In the case of a pressureless irrotational fluid, see [45,43]. The case with the null-geodesic 
equations was studied in [51]. And, the case with Boltzmann equation was considered in [52]. 

Compared with the previous works, perhaps we could emphasize the following as the new points in our work: 

(i) We present the complete sets of perturbed equations in the gauge-ready form, so that we could easily apply the 
equations to any gauge conditions which make the mathematical analyses of given problems simplest. 

(ii) We consider the most general Friedmann background with K and A. Previous studies have considered the flat 
Friedmann background only. 

(iii) We consider the most general imperfect fluid situation which includes multiple imperfect fluids with general 
interactions among them. In addition we also include minimally coupled scalar fields, a class of generalized gravity 
theories, the electromagnetic fields, the null geodesic, and the relativistic Boltzmann equation. 

(iv) In §VII A we present closed form equations which are simlar to the known ones in linear theory. 

(v) In §VII C we show that up to the second-order in perturbations the relativistic pressureless fluid coincides exactly 
with the Newtonian one. We note that suitable choices and combinations of different gauges (thus gauge-invariant 
combinations) are important to show the equivalence. 

(vi) In §VIID we have derived the large-scale (super-sound-horizon) conserved quantity to the second-order, <p v , 
directly from the differential equation governing its evolution. This consreved variable was first studied by Salopek 
and Bond 

Our equations are suitable to handle the nonlinear evolutions in the perturbative manner. If we have the solutions 
to the linear-order (see §VIIB for some examples), the evolution of second-order perturbations can be derived using 
the quadratic combination of the linear variables as sources; our basic sets of equations in §V C and some closed forms 
in §VIIA are presented with such a purpose. As long as we take such a perturbative approach our formulation in 
this work can be trivially extended to any higher-order perturbation; except for the fact that, of course, the needed 
algebra would be quite demanding. We also have shown in §VI that the gauge issue can be similarly handled even in 
such higher-order perturbation. 

Our formulation can be applied using several different methods in the following. 

(i) Quasilinear analyses using Fourier analyses as often used in the Newtonian case [44]. In this approach the 
quadratic combination of the linear-order terms will lead to the mode-mode coupling among different scales, as well 
as among different types of perturbations. 

(ii) Nonlinear backreaction. In our approach we have assumed the presence of "fictitious" background metric 
which is spatially homogeneous and isotropic. As the basic equations of Einstein gravity are nonlinear the nonlinear 
fluctuations in the metric and the matter can affect the background world model. One anticipates to recover the 
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background Fricdmann world model through averaging the more realistic lumpy world model and finding best fit to 
the idealized world model [53]. 

(iii) Fitting and averaging. Our basic equations in §IV are presented without separating the background order 
quantities from the perturbed-order ones. Thus, the equations are suitable for the operation of averaging. Using 
our formulation we could apply and check the various different averaging prescriptions suggested in the literature 
[53,54,45]. 

Our perturbative formulation would be a useful complement to the following formulations aiming to investigate the 
nonlinear evolutions of the cosmological structures. 

(i) The large-scale (long wavelength) approximation or the spatial gradient expansion studied by Salopek, Tomita, 
Deruelle and others in [46,55]. 

(ii) Cosmological post-Newtonian formulation studied by Futamase, Tomita and others in [56]. 

(iii) Rclativistic Zel'dovich approximation studied in [45,43,57]. 

(iv) General (spatially inhomogeneous and anisotropic) solutions near singularity where the large-scale conditions 
are well met; in such a situation it was shown that the spatially different points decouple and evolve separately. These 
were studied by Belinsky, Lifshitz, Khalatnikov, and others in [58]. 

Our general formulation can be used to study the following situations anticipated during the evolution of our 
universe: 

(i) We can check the limit of linear theory. Current cosmological observations can be successfully explained within 
the current standard theoretical paradigm. In that paradigm the linear perturbation theory plays significant roles 
to explain the quantum generation stage in the early universe and the classical evolution processes in the large-scale 
and in the early era. The linear theory provides a self-consistent explanation of some important aspects of the origin 
and evolution of large-scale structures. However, the limit of the linear theory cannot be estimated within the linear 
theory We expect the second-order perturbation theory could provide a meaningful ways to investigate such limits. 

(ii) We can investigate the quasilinear process in the relativistic context. In the literature it is commonly assumed 
that the relativistic linear perturbation theory is sufficient to handle the large-scale structure, and the nonlinear 
processes occur only in the Newtonian context which are often handled by the numerical simulations. The quasilinear 
evolution would be useful to investigate the transition regions between the linear and nonlinear evolutions. Our 
perturbative approach may have its own limit, because if we find the importance of second-order contributions it may 
naturally follow that higher-order contributions would become important immediately as well. Thus, we anticipate, if 
successful, the relativistic quasilinear analyses can be developed similarly as the Newtonian cosmological quasilinear 
analyses studied in [44]. 

(iii) Fate of fluctuations in the collapsing phase, and possibly through a bounce. The fluctuations of single component 
medium and the gravitational wave are described by the second-order differential equations. In the linear stage and 
in the large-scale limit, we have general solutions in eqs. (311,312,333). In an expanding phase, the C-mode is 
relatively growing and the <i-mode is decaying, thus transient. If the initial conditions (say, generated from the 
quantum fluctuations) are imposed in the early expanding phase, the d-mode parts disappear in a few e-folding 
time of the scale factor increase, thus uninteresting. The relatively growing modes both for the scalar- and tensor- 
type perturbations are characterized by the conserved amplitude of certain gauge-invariant variable, ip v (and the 
curvature variables in other gauges) and C^p, hi the large-scale limit. In the collapsing phase, however, the roles of 
growing and decaying modes are switched. In the collapsing phase the d-mode, and the vector-mode as well, grows 
quite rapidly, see eqs. (325,339,332); our solutions in §VIIB also cover the collapsing phase by considering t — ► \t\ 
with t approaching 0, see [59]. Thus, the linear perturbations grow rapidly and inevitably reach the nonlinear stage 
[6,60]. Such growths would cause the transition of our simple (spatially homogeneous and isotropic) background world 
model to the anisotropic and inhomogeneous ones studied in [58]. Although we anticipate the perturbations would 
become quite nonlinear, we hope we could investigate the transition region based on our second-order perturbation 
formulation. One simplifying fact is that in the collapsing phase the local range covered by the dynamical time scale 
<~ H^ 1 shrinks relative to the comoving scale, thus effectively the scales we are interested in satisfy the conditions 
of large-scale limit 2 . Such large-scale conditions are well met for a given comoving scale during the early evolution 
stage (near singularity) and as the background model approaches the singularity in the collapsing phase. Investigation 
of situations in the collapsing and subsequent bouncing background is left for future study; for evolutions under the 
linear assumption, sec [59] . For the general cosmological investigation near singularity, see [58] . 

(iv) It is well known that the nonlinear effect (either in the quantum generation or in the classical evolution 



2 Professor James Bardeen and Ewan Stewart have suggested that the large-scale (long wavelength) expansion or the spatial 
gradient expansion technique [46,55] would be useful to investigate such situations. 
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processes) could lead to non-Gaussian effects in the observed quantities of the CMB anisotropics and the large-scale 
galaxy distribution and motion. Maldacena has recently investigated such an effect on the CMB based on the second- 
order perturbation theory, see [47]. The first year WMAP data shows no positive detection of non-Gaussian nature 
of the CMB sky maps under a couple of non-Gaussianity tests [61]. 

(v) Evolutions in the super-horizon scale where the scale is larger than causal domain during the dynamic time 
scale. See §VII D and VII F for the conserved quantities to the second-order which were found by Salopek and Bond 
in [46]. 

(vi) In §VIIC we have shown that to the second-order a pressureless fluid with pure scalar-type perturbation 
reproduces Newtonian result. It is likely that the relativistic effect appears in the higher-order perturbation which is 
left for future investigation. 
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